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EXECUTIVE  SUMMARY 


The  Air  Force  has  recently  iaentified  a  need  for  a  low-frequency  emergency 
communications  system  using  ground  wave  propagation  in  order  to  provide 
reliable  communications  during  times  of  electromagnetic  disturbance.  To 
select  appropriate  modulation  and  error  correction  techniques  to  meet  the 
system  reliability  requirements  it  is  important  -to  know  the  detailed 
characteristics  of  the  noise  and  interference  background.  The  GWEN  system  has 
been  tentatively  assigned  to  operate  in  the  frequency  band  from  150  kHz  to  180 
kHz.  In  this  frequency  band,  there  is  reliable  ground-wave  propagation;  at 
night,  skywave  propagation  is  usually  also  seen.  The  noise  environment  during 
both  daytime  and  nighttime  conditions  is  thought  to  be  dominated  by  impulsive 
atmospheric  noise  from  lightning  strikes.  This  noise  is  known  to  have 
definitely  non-gaussian  statistics;  in  particular,  the  number  of 
high-amplitude  events  is  significantly  greater  than  predicted  from  a  Rayleigh 
amplitude  distribution.  Although  approximate  statistics  for  this  type  of 
noise  are  known  (4,5,9,29),  extensive  observations  have  not  been  performed  for 
the  frequency  range  in  question.  Accordingly,  a  series  of  observations  was 
initiated  at  NRL's  Maryland  Point  Observatory  in  southern  Maryland.  Since  the 
atmospheric  noise  is  non-stochastic  and  is  a  function  of  time-of-day,  weather 
patterns  and  season,  a  synoptic  monitoring  program  is  needed  for  full 
parameterization.  The  monitoring  site  and  period  of  observation,  however, 
included  late  spring  and  early  summer  weather  at  a  site  with  reasonably  high 
tnunderstorm  activity,  so  should  represent  an  adequate  sample  for  worst-case 
conditions.  Several  thunderstorm  fronts  passed  over  our  monitoring  site 
during  the  course  of  the  project.  Our  experiment  recorded  41  days  of  data 
quasi-continuously  between  May  19  to  July  9,  1983,  using  several  different 
equipment  configurations,  and  provides  a  satisfactory  preliminary  noise  data 
base  upon  which  to  base  GWEN  design  parameters. 

The  noise  waveforms  from  the  receiver  were  Nyquist  sampled  and  digitally 
recorded  for  later  analysis.  The  background  noise  level  was  found,  in 
general,  to  consist  primarily  of  impulsive  noise  with  non-gaussian 
statistics.  The  mean  noise  level  was  found  to  be  much  higher  between  sunset 
and  sunrise  than  during  the  day.  A  typical  nighttime  mean  noise  level  was 
lxlO”^-*W/m^/Hz,  while  a  typical  daytime  level  was  at  least  14  db  quieter. 
However,  the  daytime  level  was  at  times  significantly  higher,  especially 
during  tne  presence  of  a  local  thunderstorm  tront.  The  highest  power  impulses 
often  occurred  during  the  afternoon  period,  between  15:00  and  21:00  EDT. 

During  local  thunderstorm  activity,  impulses  with  peak  power  of 
lxlO~it+W/m2/Hz  were  observed  frequently,  while  the  most  energetic  pulse 
detected  during  our  monitoring  period  has  a  peak  power  of 
6.8x10~l4W/  m2/Hz.  The  impulsive  noise  is  well  correlated  with  active 
local  thunderstorm  systems.  The  duration  of  a  typical  lightning  noise  event 
was  of  the  order  of  0.5  seconds. 

In  addition  to  these  preliminary  conclusions,  an  important  part  of  this 
work  was  the  provision  of  an  archival  data  base  of  our  recorded  data  for 
detailed  analysis  and  comparison  with  noise  and  modulation  models  by  MITRE 
Corp  and  other  DoD  authorized  users. 


SYNOPTIC  OBSERVATIONS  OF  THE  RADIO  NOISE  BACKGROUND 
IN  THE  FREQUENCY  RANGE  150-180  kHz  TO  PROVIDE 
DESIGN  DATA  FOR  THE  LOW  FREQUENCY 
COMMUNICATIONS  DESIGNER 

EXPECTLD  RADIO  NOISE  brt(X  GROUND  IN  THE  FREQUENCY  RANGE 
OF  150  to  180  kHz 

The  properties  of  the  radio  frequency  range  below  00  kHz  (VLF)  have  been 
extensively  studied  (23),  but  those  of  the  L F  range  (30-300  kHz)  have  been 
much  Less  so,  although  a  general  data  base  has  been  compiled  {29).  In 
general,  LF  radio  propagation  can  be  expected  to  share  many  of  the  properties 
of  the  VLF  spectrum  modified  by  differences  in  propagation  in  the  two 
frequency  domains.  Our  investigations  generally  supported  this  conclusion. 

An  estimate  of  the  magnitude  of  the  average  noise  can  be  obtained  from  the 
data  in  C.C. I.R.  Report  322  (29).  For  the  location  of  the  Maryland  Point 
Observatory  the  curves  in  this  report  give  expected  noise  levels  at  163  kHz 
corresponding  to  noise  factor  (FaJ  values  between  95  dB  for  the  daytime  low 
to  a  high  nighttime  peak  of  121  dB.  In  the  units  used  in  this  report,  this 
corresponds  to  6.bxlO“^W/m2/Hz  and  2. 6xlQ“^W/m2/hz  respectively. 

In  parameterizing  the  LF  noise  background,  there  are  two  important 
considerations;  the  source  of  the  natural  noise  background,  ana  the 
propagation  conditions  from  the  noise  source  to  the  receiver.  A  third 
important  consideration  in  this  frequency  range  is  man-made  interference,  but 
no  systematic  study  was  made  of  this  for  this  report.  Systematic  models  have 
been  made  of  the  propagation  characteristics  and  noise  background  for 
VLF(11, 13, 15) ,  but  cannot  be  extended  to  LF  without  modification. 

Atmospheric  radio  noise  in  the  VLF  band  is  primarily  a  result  of  lightning 
discharges  occurring  in  the  earth's  atmosphere  (5,6,7,17,19,21).  For 
electric-field  antennas  located  in  dry,  windy  environments  precipitation 
static  can  become  an  important  noise  contributor.  On  occasions  solar 
disturbances  can  cause  a  noise  background  increase  known  as  Sudden  Enhancement 
of  Atmospherics.  Nyquist  noise  associated  with  the  receiver,  or  with  the 
terrestrial  or  extraterrestrial  background,  is  generally  negligible,  unlike 
the  situation  at  microwave  frequencies.  Propagation  upward  or  downward 
through  the  ionosphere  has  a  high  enough  attenuation  to  make  extraterrestrial 
noise  a  negligible  factor  in  the  150-180  kHz  band,  although  signals  can  be 
propagated  through  the  earth-space  interface( 2) . 

The  large  impulsive  type  currents  associated  with  lightning  phenomena 
cause  electromagnetic  radiation  whicn  is  propagated  to  the  receiver;  this 
radiation  has  been  investigated  by  Watson  Watt  and  others  (18,20,21,22,25). 

The  initiating  predischarge  or  'leader'  stroke  develops  over  a  period  of  a  few 
milliseconds.  It  consists  of  a  series  of  discrete  sub-leaders  occurring  at  a 
rate  of  one  every  .025-. i  millisecond.  Ihe  total  time  required  for  Che  leader 
to  reach  the  ground  is  .5-1  millisecond.  Following  the  pre-discharge  leader, 
a  single  ground-to-cloud  discharge  stroke  occurs  lasting  .1  milliseconds  with 
a  peak  intensity  of  20,000  amperes.  Frequently,  'long'  discharges  can  occur 
consisting  of  several  strokes  each  lasting  .1  milliseconds,  separated  by  an 
average  period  of  40  milliseconds;  chus  a  single  stroke  can  extend  over  a  time 
interval  of  several  hundred  milliseconds.  The  total  power  in  a  single 
discharge  lasting  100  milliseconds  is  500-1000  megawatts,  much  of  which  goes 
into  ionizing  the  atmosphere  and  generating  acoustic  noise.  It  has  been 
estimated  that  as  much  as  1  megawatt  remains  for  generating  the  RF  field 
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(6,23).  The  radio  frequency  noise  Chat  is  generated  is  impulsive  in 
character  (3),  with  a  duration  that  corresponds  to  the  physical  duration  of 
the  lightning  stroke.  The  peak  frequency  of  the  radio  frequency  noise  is  in 
the  broad  range  of  5-50  kHz;  estimated  signal  strengths,  normalized  to  a 
distance  of  1  km.,  can  be  as  high  as  2xl0-7W/nr /Hz  at  a  carrier  frequency 
of  5  kHz.  The  appropriate  estimate  for  the  vicinity  of  100  kHz  is 
2x10-12^/3,2 aC  a  distance  of  1  km,  and  is  estimated  to  gradually 
decrease  with  increasing  frequency  at  a  rate  of  3  dB/octave.  This  is  still  a 
highly  significant  source  of  noise. 

The  occurrence  of  lightning  discharges  varies  from  place  to  place  on  the 
earth(lO).  Lightning  discharges  are  very  rare  in  Arctic  and  Antarctic  regions 
and  can  be  as  low  as  1  discharge  per  second  in  temperate  regions  or  as  high  as 
100  per  second  in  the  tropics.  Ac  any  one  time  there  are  approximately  2000 
thunderstorms  in  progress  worldwide  (2).  In  the  temperate  regions  especially, 
the  thunderstorm  incidence  is  a  strong  function  of  the  season  of  the  year. 

Ground  wave  propagation  and  sky  wave  propagation  (ionospheric  bounce)  are 
the  principal  modes  of  propagation  that  take  place  at  LF  frequencies.  The 
ground  wave,  which  is  the  basis  of  the  GWEN  system,  is  a  very  effective 
propagation  mode  at  165  kHz  due  to  refraction  at  the  earth's  surface;  for 
example,  at  a  distance  of  1000  km  the  ground  wave  should  be  attenuated  by  only 
about  20  dB  more  than  the  l/r2  geometrical  attenuation.  Ionospheric 
propagation  is  considerably  more  variable;  however,  a  clear  day-night  effect 
is  usually  present,  with  the  sky-wave  propagation  being  much  more  efficient 
from  sunset  to  sunrise  than  during  the  daytime.  This  is  due  to  a  difference 
in  ionospheric  reflection  coefficient  Chat  is  somewhat  complex  to  predict,  but 
is  due  basically  to  the  reflection  occurring  at  lower  altitudes  during  the  day 
where  the  higher  density  causes  a  decrease  in  effective  ionospheric 
conductivity.  The  net  power  delivered  to  the  receiver  by  skywave  propagation 
from  distant  thunderstorms  can  be  appreciable  if  there  is  a  large  thunderstorm 
front  being  propagated.  In  addition,  a  focusing  effect  is  possible  giving 
path  losses  significantly  less  than  l/r2.  Waveguide -mode  propagation, 
although  important  at  VLF(l),  is  not  expected  to  be  a  significant  effect  in 
the  150-180  kHz  range. 

A  typical  ground  station  is  thus  expected  to  observe  two  noise  components; 
local  noise  from  nearby  thunderstorms  within  500-1000  miles,  propagated  via 
the  ground  wave,  and  a  component  due  to  distant  thunderstorms  propagated  via 
sky  wave.  Since  the  local  component  consists  of  contributions  from 
appreciably  fewer  lightning  strokes,  it  might  be  expected  that  it  would  be  of 
significantly  more  impulsive  character  than  the  sky-wave  component  which  might 
consist  of  many  pulses  overlapping  in  time  to  create  a  quasi-continuous 
background.  Improved  sky-wave  propagation  will  cause  the  background  level  to 
increase  as  progressively  more  distant  thunderstorms  are  able  to  contribute  to 
the  received  RF  field.  The  impulsive  component,  on  the  other  hand,  should  be 
much  more  sensitive  to  local  weatner  where  the  RF  field  is  dominated  by 
individual  lightning  strokes. 

The  robustness  of  a  receiver  system  in  the  presence  of  background  noise 
depends  on  both  the  characteristics  of  the  noise  background  and  on  the 
modulation  characteristics  of  the  desired  signal;  our  present  work  did  not 
investigate  this  subject  in  depth.  Work  on  receiver  noise  sensitivity  has 
been  performed  by  Watt  et  al  and  others  (14,24);  there  is  a  fairly  extensive 
literature  on  characteristics  of  non-gaussian  noise  expected  at  these 
frequencies  (16,26,27,28).  The  present  investigation  focused  on  obtaining  a 
valid  statistical  data  base  for  further  analysis,  and  on  determining  basic 
noise  parameters  in  the  frequency  range  considered. 


EXPERIMENTAL  PROCEDURE 


To  study  atmospheric  noise  realistically  requires  a  high  dynamic-range 
receiver  that  can  produce  wide-bandwidth  samples  of  natural-background  noise 
over  a  selected  period  of  observation.  In  addition,  the  need  for  rapid 
results  required  the  use  of  equipment  available  "off-the-shelf".  From 
available  equipment,  it  was  possible  to  make  a  working  system  which  consisted 
of  a  magnetic  pickup  loop,  connected  through  preamplifiers  to  commercial 
(Hewlett-Packard)  wave  analyzers,  which  were  set  to  act  as  a  single-sideband 
superheterodyne  receiver  with  adjustable  bandwidth  and  center  frequency.  The 
bandwidth  was  normally  set  to  3.1  kHz,  the  largest  value  available. 

Undetected  i.f.  noise  instead  of  envelope  data  was  recorded  to  maintain  full 
flexibility  in  later  coherent  signal  analysis.  Envelope  and  other  information 
was  created  by  computer  post-processing.  During  the  latter  portion  of  the 
experiment,  a  computer-controlled  attenuator  was  used  to  set  the  system  gain 
in  steps  of  approximately  20  dB  to  extend  the  system  dynamic  range.  The 
atmospheric  noise  level  during  the  quietest  daytime  conditions  was  about  5  dB 
higher  than  the  output  of  the  first  preamplifier  at  the  loop  site.  The  video 
output  of  the  wave  analyzer  was  digitally  sampled  at  a  rapid  rate  and  recorded 
on  a  computer  tape.  The  recording  system  was  controlled  by  a  Hewlett-Packard 
2116C  computer  system  that  had  previously  been  used  in  several  similar 
real-time  equipment  control  projects.  The  block  diagrams  of  the  equipment 
are  shown  in  Figs.  1-4,  and  the  data  control  block  diagram  is  shown  in 
Fig.  5.  Two  channels  of  data  were  recorded  with  12-bit  precision;  various 
configurations  were  used  throughout  the  course  of  the  experiment  as  detailed 
in  Figs.  1-4  to  record  different  signal  levels  to  enhance  the  dynamic  range. 
The  data  log  is  transcribed  in  Tables  I-III.  A  dynamic  range  of  120  dB  was 
desirable;  the  12-bit  analog-to-digital  converters  used  had  an  effective 
dynamic  range  without  changing  settings  of  about  60  dB,  as  did  each  wave 
analyzer  used.  The  second  channel  was  usually  set  at  about  40  dB  lower  gain 
to  give  a  total  range  of  100  dB.  The  actual  data-taking  format  was  controlled 
in  a  versatile  manner  by  a  real-time  HP2116  computer  programmed  in  Forth.  A 
block  diagram  of  the  data-taking  program  is  shown  in  Fig.  5;  listings  of  the 
code  are  in  Appendix  B,  but  it  should  be  noted  that  a  full  explanation  of  the 
code  is  not  attached. 

The  data-taking  routine  operated  approximately  as  follows;  at  start  of 
execution  of  the  Forth  word,  values  for  the  various  parameters  detailed  in  the 
identification  section  were  typed  in  from  the  system  terminal.  These  included 
frequency,  bandwidth,  location,  sample  interval,  etc.  The  start  time  was 
obtained  from  the  U.S.  Naval  Observatory  telephone  time  service  and  was  used 
to  start  the  internal  time-keeping  of  the  computer,  which  consisted  of  a 
crystal-controlled  0.1  second  interrupt  generator,  with  an  epoch  accuracy  of 
about— 0.2  seconds.  When  time  for  a  sample  to  be  taken  was  at  hand,  control 
was  transferred  to  a  sampling  routine.  This  routine  recorded  data  at  a  rate 
which  was  controlled  by  interrupts  generated  by  an  external  synthesizer,  which 
was  set  at  approximately  the  Nyquist  rate  for  the  bandwidth  used.  For  the 
6  kHz  rate  used  for  most  of  the  sampling,  the  sample  interval  was 
167  microseconds.  The  sample  from  the  second  channel  was  recorded  close  to 
simultaneously  with  that  from  the  first  channel;  there  was,  however,  a 
23  microsecond  delay  due  to  hardware  constraints  in  the  multiplexer.  A  tape 
record  consisted  of  5000  samples  from  each  channel,  covering  a  total  interval 
of  0.833  seconds,  together  with  identification  information.  Some  of  the  early 


data  was  recorded  wich  a  synthesizer  frequency  of  10  kHz,  resulting  in  a 
sample  interval  of  100  microseconds  and  a  record  length  of  0.5  seconds.  Later 
data  was  recorded  in  groups  of  ten  records  to  ensure  good  sampling  of  the 
medium-term  statistics  of  the  lightning-associated  noise.  There  is  a  data  gap 
of  0.022  seconds  it 0.001  seconds  between  tape  records  within  a  group  to  allow 
for  program  overhead  in  recycling  for  the  next  record.  A  major  limitation  in 
the  amount  of  data  recorded  was  the  volume  and  expense  of  the  magnetic  tapes. 

A  limit  was  somewhat  arbitrarily  set  of  recording  one  magnetic  tape  per  day; 
at  our  800  bpi  recording  density  this  meant  that  we  could  record  about  1000 
seconds  per  day  total,  or  about  li  of  the  time. 

For  the  majority  of  the  data,  both  channels  were  recorded  at  the  same 
nominal  center  frequency,  with  all  parameters  the  same  except  for  the  gain 
setting.  The  question  then  arises  to  what  extent  the  two  channels  can  be 
combined  to  equivalence  a  single  channel  with  larger  dynamic  range.  The 
answer  for  our  experiment  is  that  for  purposes  of  determining  the  noise 
envelope,  data  from  both  channels  can  be  combined  freely.  However,  both 
channels  cannot  be  combined  in  a  fully  coherent  fashion  primarily  because  the 
bandpass  characteristics  and  center  frequencies  were  not  identical. 

Calibration  of  the  URM-6  loop  gain  was  performed  by  the  method  described 
in  Jean  et  al  (12)  to  establish  the  relationship  between  the  signal  input  to 
the  test  port  of  the  loop  and  field  strength.  An  overall  system  calibration 
was  performed  by  introducing  a  known  signal  level  into  the  test  port  of  the 
loop  and  recording  the  output  for  the  data-taking  system.  Because  the  antenna 
and  the  preamplifier  were  operated  as  an  integral  unit,  no  sensitivity  figure 
was  determined  for  the  antenna  alone;  but  the  overall  system  calibration  was 
determined  to  be:  1  unit  of  power  *  2. 78X10“21w/m2/nz>  where  one  unit  of 
power  is  defined  as:  one  count2  on  channel  one,  with  the  digital  attenuator 
set  at  a  value  of  1  and  all  system  gains  set  as  indicated  in  Fig.  3. 

It  should  be  noted  that  all  signal  strengths  in  this  report  are  given  in 
units  of  power  spectral  density,  watts  per  square  meter  per  unit  bandwidth. 
This  type  of  measure  is  appropriate  for  impulse-type,  or  wideband,  noise.  The 
power  units  are  easily  convertible  to  field  strength  in  microvolts  per  meter 
by  use  of  the  formula: 

9 

P  ■  E“ 

R 

where  P  *  Power  density 
E  *  field  strength 

R  *  the  impedance  of  free  space,  376.3  ohms 

The  overall  frequency  response  of  the  system  was  calibrated  by  injecting  a 
signal  into  the  antenna  and  measuring  the  level  at  the  input  to  the 
analog-to-digital  converter.  The  result,  which  is  the  overall  system 
response,  is  shown  in  Figure  7.  Figure  7  is  for  the  HP3590A  wave  analyzer, 
which  was  normally  considered  our  primary  channel.  At  a  nominal  observing 
frequency  of  165  kHz,  a  frequency  of  163.2  kHz  is  converted  to  zero  frequency, 
and  the  mean  observed  frequency  is  about  164.8  kHz.  A  response  curve  for  the 
wave  analyzer  alone  is  shown  in  Figure  6. 


OBSERVATIONS 


Because  of  the  rapid  time  scale  of  the  project,  the  final  equipment 
configuration  was  decidea  upon  while  the  experiment  was  in  progress.  Four 
equipment  configurations  were  used  for  appreciable  periods  of  time.  These  are 
labeled  Systems  A,  B,  C  and  D  and  are  set  forth  in  block  diagram  form  in  Figs. 
1-4.  While  the  latter  configurations  were  in  general  an  improvement  over  the 
earlier  ones,  the  earlier  data  was  judged  worth  retaining.  The  data  from  the 
period  before  13  June  (Systems  A  and  B)  is  most  useful  for  purposes  of 
determining  the  average  noise  statistics,  that  during  the  period  13-27  June 
(System  C)  for  determining  peak  statistics,  while  that  taken  on  28  June  and 
after  (System  D)  should  be  useful  for  both  purposes. 

The  experiment  was  initially  set  up  in  Building  209  at  the  Naval  Research 
Laboratory,  Washington,  D.C.  However,  it  was  quickly  determined  that  various 
sources  of  man-made  local  interference,  primarily  from  computers,  were 
dominating  the  output,  so  the  base  of  operations  was  transferred  to  our  radio 
astronomy  field  station  at  Nanjemoy,  Maryland.  This  station  is  located  on  the 
Potomac  River  in  a  primarily  rural  area  about  40  miles  south  of  Washington. 

Its  approximate  geographic  latitude  is  38  22. 4'  North,  and  its  longitude  is 

77  13.9'  West.  Local  interference  was  not  a  significant  problem  at  this 

site.  There  were,  as  expected,  several  interfering  carriers  present  in  the 
general  frequency  band;  the  primary  observation  frequency  of  165  kHz  was 
chosen  to  minimize  this  problem  after  monitoring  of  the  frequency  band. 
Observations,  undertaken  quasi-continuously  during  the  period  13  May  -  8  July 
1983,  were  recorded  on  a  series  of  magnetic  tapes.  The  background  was 
monitored  at  one-to-ten  minute  intervals  during  almost  all  of  this  interval; 
occasional  exceptions  were  due  to  configuration  changes,  equipment  failure  or, 
most  often,  power  failure  at  the  site  due  to  thunderstorm  activity.  During 
this  time,  there  was  a  fairly  typical  assortment  of  Washington-area 
early-summer  weather,  including  several  periods  of  local  thunderstorm  activity 
as  well  as  fair-weather  periods.  The  weather  environment  was  monitored  in  a 
quantitative  manner  by  tne  installation  of  a  rain  gauge  at  the  site,  as  well 
as  by  daily  videotape  recording  of  the  national  weather  pattern,  including  the 
digital  radar  map  that  is  an  excellent  indicator  of  precipitation  activity. 

In  analyzing  data  from  a  specific  period  from  this  project,  it  is  vital 
that  the  user  refer  to  the  data-related  tables  at  the  back  of  this  report,  the 
synoptic  atlas  of  LF  noise  levels.  There  are  three  data-related  tables. 

Table  I,  entitled  List  of  Significant  Changes  and  Settings  in  Data-taking 
Procedure,  lists  all  the  generic  changes  in  data-taking  procedure  and  is 
important  to  consult  to  understand  the  data  recorded.  Table  II,  entitled  List 
of  Unique  Characteristics  of  Individual  Tapes ,  records  non-standard  features 
of  particular  tapes.  Table  III,  the  Data  Tape  Log,  includes  start  time, 
interval,  and  all  otner  data  parameters  that  were  changed  on  a  daily  basis,  as 
well  as  a  summary  of  the  weather  information. 

A  significant  purpose  of  this  experiment  was  to  create  an  archival  data 
base  of  LF  noise  data  that  could  be  used  for  signal  processing  simulations 
during  the  GWEN  design  effort.  NRL  is  therefore  willing  to  supply  copies  of 
this  tape  data  base  at  nominal  duplication  cost  to  autnorized  groups.  For 
further  information  contact: 

S.H.  Knowles 
Code  4 183 A 

U.b.  Naval  Research  Laboratory 
Washington,  D.C.  20375 
ph.  202-767-3010 


theory  of  analysis  of  undetected  wideeand  samples  of  lf  noise 


To  maintain  flexibility  in  later  analysis,  date  was  recorded  as  a  wideband 

undetected  data  stream.  Since  data  of  this  type  is  not  in  wide  use  in  the 

field  of  LF  signal  analysis,  it  is  appropriate  to  discuss  some  basic  precepts 
for  analysis.  A  6  khz  sampling  rate  is  assumed. 

The  detected  voltage  at  165  kHz  in  a  3100  Hz  bandpass  is  down-converted  to 
a  base  band  at  15o0  Hz  center  frequency  and  then  Nyquist  sampled  at  6000  Hz 
for  a  period  of  833  milliseconds.  The  instantaneous  voltage  can  be  described 
as 

vlt)  =  V(t)  3 in  (2»  f  t  +  <»(t)  ) 

where  v(t)  is  the  slowly  varying  voltage  of  the  noise  envelope,  f  is  the 
frequency  of  the  center  of  the  observing  band  (1550  Hz)  and  <t>  represents  the 

phase  of  the  received  signal.  The  average  received  noise  power,  P,  during  a 

sample  period  of  T  seconds  is 

P  ®  (  V( t )2  sin2  (  2  r  f  t  +  <P  (t)  )  y 

Since  the  envelope  power  is  a  slowly  varying  function  of  time  and  is 
physically  uncorrelated  with  the  band  center  frequency  f, 

P  -  <V(t)2)  ,<sin2  (2»f  t  +b(t)  )> 

which  simplifies  to 

P  -  JL_  <  V(t)2  > 

2 

where 

<  V(t)2  >  *■  .1 
T 

We  see  that  the  mean  observed  power  P  is  equivalent  to  the  average  envelope 
power  V(t)2. 

In  addition  to  the  average  power  in  each  833  millisecond  sample  it  is  also 
possible  to  calculate  the  voltage  deviation  Vd  which  is  related  to  the 
impulsiveness  of  the  incoming  signal.  It  is  defined  in  terms  of  the 
instantaneous  envelope  voltage  as  follows, 

Vd  =  20  Log 

where 

<  V'(tJ2>  =  l  f  v(t)2  dt 

T  Jo 

<(V'(t)  y  *  if  lv(t)l  dt 


The  quantity  <V'lt;2>  can  be  easily  derived  from  our  expression  for  the 
detected  mean  power  P 


p  <V(t)2>  =  <V'(t)2> 

2 

The  quantity  <)V'(tJl>  is  obtained  from  our  measured  quantity  v(t)  by 
the  relationship 

OV'itJl)*  =  <  IVUJ  sin  (2ir  f  t  +  <S>  It))!)* 

*  <jV(tJl)>  <(lsin(2rrf  t  +  <J>  (t))l)‘ 

=  2/t  *  <IV(t)l> 


AMPLITUDE  AND  TIME  PROBABILITY  DISTRIBUTION 


The  Amplitude  Probability  Distribution  (APD.)  represents  the  percentage  of 
time  that  the  noise  power  exceeds  a  given  threshold  level.  If  the 
instantaneous  noise  power  is  defined  as  P(t)  =  Vit)  ,  we  can  construct  the 
function  p(t)  such  that 


p(t) 


1  for  P(U  >  P0 

0  for  P(t; <  P0 


where  PQ  represents  the  threshold  level.  We  then  integrate  this  function 
over  the  sample  interval  to  determine  the  percentage  of  time  that  the  signal 
exceeds  the  value  P 


APD(P0) 


p(t)  dt 


A  second  measure  of  the  observed  noise  is  the  Time  Probability 
Distribution  1TPD).  This  function  can  be  obtained  from  the  data  in  the 
following  manner.  We  select  a  threshold  power  level  P0  and  from  the  arrival 
times  of  the  noise  impulses,  we  construct  the  function 


P(t)  *  1  for  P(t)  >  PD 

*  0  for  P(t)  <  PQ 


This  eliminates  all  information  regarding  the  power  in  each  impulsive  event 
whose  distribution  has  already  been  defined  by  the  APD.  From  the  function 
P(t)  we  compute  the  normalized  autocorrelation  function 


A(r) 


I 


L 

r. 


P(t).P(t-r).dt 


Pit! . Pit) .dt 


This  function  has  the  desired  property  that  it  counts  up  the  number  of 
occurrences  for  which  the  pulse  spacings  above  the  threshold  power  P0  are 
separated  by  a  time  interval  T.  By  normalizing  the  function  Ad;  computed  for 
a  particular  threshold  by  the  total  number  of  crossings  of  the  signal  at  P0> 
the  relative  frequency  of  occurrences  of  pulses  at  a  separation  T  can  be 
determined.  The  TPU  may  be  derived  from  AiTl  by  integrating  A(T) 
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t 


A(t  Me 


TPD(T  >  t)  ■  100&  Jl 


Jt 


L 


A(tMc 


0 


which  now  gives  the  percentage  of  time  that  the  pulse  spacings  at  a  given 
threshold  exceed  the  time  t. 


INTERPRETATION 


The  average  background  power  at  165  kHz  on  June  8-9,  1983  is  displayed  in 
Figure  8.  The  background  is  highly  non-uniform  over  a  period  of  24  hours  and 
also  varies  from  day  to  day.  Sunrise  and  sunset  occurred  at  approximately 
6:00  and  20:30  EST;  and  as  Figure  8  indicates,  these  diurnal  events 
consistently  correspond  to  rapid  changes  in  the  average  noise  background.  The 
nighttime  level  was  usually  14  dB  higher  than  the  'sytime  level. 

An  analysis  of  the  amplitude  and  time  probability  distributions  provides  a 
valuable  tool  for  studying  the  statistical  characteristics  of  the  background 
noise.  For  example,  if  the  noise  has  a  Rayleigh  distribution  we  should  see  an 
APD  with  the  form 


N(P>PS)  *  Constant  x  exp(-.694  P/Ps) 

where  P  represents  the  threshold  power  level  and  Ps  is  the  half-power  width 
of  the  distribution  function.  We  have  selected  three  typical  samples  of  data 
to  represent  daytime,  nighttime,  and  thunderstorm  noise  properties  (Figure  9) 
to  be  used  in  our  analysis.  A  comparison  of  the  theoretical  curve  with  a 
typical,  daytime  APD  from  a  single  .833  second  sample  (Figure  10)  shows  that 
the  noise  follows  a  Rayleigh  distribution  below  a  threshold  of 
1 . lxl0"^^W/m^Hz ,  a  value  that  is  6  times  the  average  power  level  for  the 
sample.  Above  this  threshold,  there  are  noticeably  more  impulsive  events  than 
predicted  from  the  Rayleigh  statistics.  In  Figure  10  we  also  show  the  APD  for 
a  .833  second  sample  for  the  nighttime  period  on  the  same  day.  Here  we  see 
tnat  the  data  is  Rayleigh  distributed  below  a  thresnold  of  lxl0“^-7W/m2/Hz, 
once  again  a  value  that  is  6  times  the  average  power  for  that  sample,  and  is 
noticeably  non-Rayleigh  above  this  threshold.  The  average  noise  level  differs 
by  a  factor  of  5-10  between  these  two  samples,  which  accounts  reasonably  well 
for  the  difference  between  the  two  threshold  levels.  A  comparison  of  the 
day/night  APDs  shows  that  the  nighttime  distribution  is  significantly 
flatter.  The  number  of  events  exceeding  6  times  the  average  noise  level  is 
much  higher  in  the  nighttime  sample  than  in  the  daytime  sample. 

A  comparison  of  the  day/night  APD  with  the  APD  from  a  sample  taken  during 
a  thunderstorm  on  May  26  (Fig.  10)  shows  a  much  different  character,  however. 
The  thunderstorm  APD  shows  a  greater  frequency  of  impulsive  events  between  2 
and  17  times  the  average  background  level  than  expected  from  Rayleigh 
statistics . 

The  predicted  TPD,  assuming  that  the  arrival  rate  is  completely 
uncorrelated,  can  be  expressed  as  a  gaussian  process 

P( t  >  T)  =  exp (-NT) 

where  N  is  the  average  arrival  frequency.  We  may  interpret  the  observed  TPDs 
by  comparing  them  to  this  model  of  a  completely  uncorrelated  pulse  arrival 
distribution  which  yields  a  straight  line  in  Figure  11  with  a  slope  of  -1.  We 
see  tnat  when  the  threshold  is  set  to  zero  as  in  curve  'A1,  the  arrival  rates 
remain  uncorrelated  for  time  intervals  shorter  than  400  ms.  Since  the  full 
record  length  is  only  833  ms,  we  are  clearly  unoersampling  pulse  intervals 
separated  by  times  greater  than  400  ms.  By  increasing  the  record  length  to 
8.3  seconds  this  limitation  can  probably  be  avoided. 


We  have  computed  the  TPD  for  two  additional  thresholds  (4  and  C)  which 
correspond  to  levels  at  3.7  and  1.4  dB  below  the  highest  detected  peak  in  each 
sample.  As  we  see  from  the  daytime  and  nighttime  distributions,  as  the 
threshold  is  increased,  the  TPD  deviates  markedly  from  the  case  of 
uncorrelated  behavior  represented  by  the  dashed  lines  with  slopes  of  -1.  This 
deviation  occurs  over  an  interval  of  8  ms  to  loO  ms.  In  particular,  we  see 
that  at  the  highest  threshold,  78%  of  all  observed  pulse  spacings  ought  to  be 
separated  by  80  ms  whereas  for  the  daytime  sample  (curve  C)  over  98%  of  the 
pulse  spacings  exceed  this  interval  implying  a  high  degree  of  correlation.  As 
we  see  from  the  data  itself  (Figure  9J,  the  TPD  is  biased  since  the  number  of 
threshold  crossings  is  only  5  as  compared  to  36  for  the  3.7  dB  level  and  3000 
for  the  zero  level  thresholds. 

For  threshold  ri,  the  TPD  shows  that  impulsive  events  with  intervals  less 
than  16  ms  are  more  strongly  correlated  in  their  arrival.  There  are  more  of 
these  events  observed  than  expected  from  a  random  distribution.  Events  having 
spacings  longer  than  this  are  very  nearly  completely  uncorrelated  over  the 
range  from  16ms  to  300  ms. 

For  the  nighttime  sample,  the  distribution  is  markedly  different.  Once 
again,  for  intervals  greater  than  400  ms  the  undersampling  of  tne  data  becomes 
a  problem.  Generally,  the  distribution  between  16ms  and  30ms  is  uncorrelated 
at  the  two  highest  thresholds  as  was  the  case  for  the  daytime  sample.  For 
intervals  less  than  16ms,  we  once  again  notice  a  departure  from  completely 
uncorrelated  behavior,  this  time  in  the  direction  of  fewer  such  events  than 
predicted. 

Evidently,  the  daytime  and  nighttime  noise,  based  on  a  single  .833  second 
sample  of  average  conditions  during  each  period,  is  quite  similar.  The  pulse 
interval  distributions  tend  to  be  uncorrelated  between  16  and  300ms.  For 
times  shorter  than  this,  the  behavior  tends  to  be  non-random  and  correlated. 
This  can  be  understood  if  the  typical  duration  of  a  single  impulsive  event 
were  ^16ms  so  that  we  are  resolving  an  individual  burst  of  noise.  An 
examination  of  Figure  9  shows  that  this  is  indeed  the  case. 

In  Figure  10  we  have  also  determined  the  TPD  for  a  representative  sample 
of  data  taken  during  a  thunderstorm  on  May  26.  We  see  that  for  spacings 
shorter  than  110  ms,  the  noise  impulses  appear  to  be  uncorrelated  at  all  three 
thresholds,  however  for  intervals  between  110ms  and  300ms,  the  distribution  is 
noticeably  correlated.  An  inspection  of  Figure  9  shows  that  the  noise  tends 
to  be  clumped  at  intervals  spanning  the  range  110  to  300  ms.  This  behavior  is 
far  more  pronounced  than  the  daytime  or  nighttime  events.  This  distribution 
can  be  explained  by  lightning  discharges  consisting  of  multiple  strokes  at 
intervals  of  .1  to  .3  second,  each  consisting  of  a  random  process,  perhaps  the 
propagation  of  pre-discharge  leaders.  Once  again,  it  should  be  kept  in  mind 
that  the  above  discussion  represents  the  analysis  of  a  single  .833  second  data 
sample  for  a  lightning  discharge.  To  what  extent  this  is  typical  of  all  such 
discharges  will  need  to  be  determined  by  studying  a  larger  sample  of 
thunderstorms  during  their  entire  liretimes. 

The  majority  of  the  data  recorded  and  analyzed  was  at  a  nominal  frequency 
of  165  kHz.  However,  between  June  13  and  June  27  three  days  of  data  were 
recorded  with  the  second  channel  at  150  kHz  and  seven  days  at  180  kHz.  An 
examination  of  this  data  shows  no  significant  difference  between  the  two 
channels  indicating  no  dependence  of  the  noise  characteristics  on  frequency 
between  150  and  180  khz. 


«•. 
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TaBLE  1 


LIST  OF 
5/10/83 
5/13/83 
5/19/83 


5/23/83 


5/26/83 


GWEN  PROJECT 


SIGNIFICANT  CHANGES  AND  SETTINGS  IN  DATA-TAKING  PROCEDURE 

-  Very  preliminary  tape  -  not  kept  as  part  of  data  base 

-  One-channel  test  tape. 

-  Two-channel  setup  with  one  wave  analyzer  (H-P  3590A) ,  two  operational 
amplifiers  for  different  sensitivities.  Gain  settings  as  follows: 

H-P  46 1A  preamp  -  40  dB  gain 

H-P  3590A  wave  analyzer:  input  voltage  -  0.03  volts 

range  -  -50  dB. 

center  frequency  -  165.00  kHz 
bandwith  -  3100  Hz 
mode  -  upper  side  band 
strip  chart  -  5  volts  full  scale 
H-P  400E  power  meter  -  1  volt 
H-P  606B  cal.  generator  -  -120  dBm  output 

operational  amplifier  gains:  channel  #0  -  X2 

channel  #1  -  X20 

-  New  standard  gain  settings  to  prevent  saturation: 

H-P  461a  preamp  -  20  dB  gain 

H-P  3590A  wave  analyzer:  input  voltage  -  0.3  volts 

range  -  -50  dB. 

center  frequency,  bandwidth  and 

sideband  unchanged 

-  Equipment  configuration  changed  to  include  two  wave  analyzers, 
to  increase  dynamic  range.  Sample  rate  decreased  from  10  kHz  to 
6  kHz,  increasing  length  of  one  record  to  0.83333. . .seconds. 

Monitor  routine  added  to  ensure  that  possible  "hang-up"  of 
analog-to-digital  converter  is  monitored. 

New  standard  gain  settings  as  follows: 

H-P  461A  preamp  -  40  dB  gain 

input  voltage  -  0.1  volts 
range  -  -50  dB. 

center  frequency  -  165.00  kHz 
bandwith  -  3100  Hz 
mode  -  upper  side  band 
A/D  assignment  -  channel  J/0 
input  voltage  -  10  volts 
range  -  -50  dB. 

center  frequency  -  165.00  kHz 
bandwith  -  3000  Hz 
mode  -  lower  side  band 
A/D  assignment  -  channel  ;/i 
strip  chart  -  5  volts  full  scale 
H-P  400E  power  meter  -  1  volt 
H-P  606B  cal.  generator  -  -120  dBm  output 


H-P  3590A  wave  analyzer: 
(high-gain  channel) 


H-P  310A  wave  analyzer: 
(low-gain  channel) 
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operational  amplifier  gains:  channel  #0  -  X2 

channel  ?/l  -  X5 

note:  input  voltage  ranges  of  one  or  both  wave  analyzers  may  be 

lower  by  a  range  (X3)  in  first  day  or  so. 

5/27/83  -  New  calibration  routine  :  rec.  #1  -  -120  dBm 

rec.  It2  -  -90  dBm 

5/31/83  -  Software  changed  to  sample  a  group  of  10  records  at  a  time  in 
order  to  better  sample  longer  r  s.  Total  group  length  is  now 
increased  to  8. 33333. .. seconds ,  with  about  0.022  sec. 
inter-record  gap,  ano  inter-group  interval  increased  to  10  min. 
from  1  min. 

6/1/83  -  Standard  inter-group  interval  increased  to  12  minutes. 

6/9/83  -  Removed  230  kHz  low-pass  filter  and  changed  from  AIL  to  H-P 

synthesizer. 

6/13/83  -  Digitally-controlled  attenuator  installed  to  control  input 

range  automatically.  Most  gain  settings  remain  unchanged.  Both 
wave  analyzers  are  now  set  to  input  voltages  of  0.1  volt  and 
ranges  of  50  dB.  The  310A  (channel  #1)  is  now  set  at  an 
alternate  frequency. 

Constants  for  attenuator  algorithm: 

upscale  trigger  -  1024  counts 

downscale  trigger  -  100  counts 

wait  for  downscale  -  10000  samples 
H-P  310A  center  frequency  set  to  150  kHz.  No  calibration  records 
recorded  for  remainder  of  data. 

6/14/83  -  Calibration  tape  recorded  to  calibrate  impulse  response  of  wave 
analyzers,  and  to  measure  inter-record  interval. 

6/17/83  -  H-P  310A  center  frequency  changed  to  180  kHz. 

6/23/83  -  Recorded  calibration  tape  for  absolute  amplitude  calibration. 

6/28/83  -  Algorithm  on  switched  attenuator  changed  to  prevent 

over-reduction  of  range,  and  ensure  wide  dynamic  range.  In  the 
new  configuration,  both  wave  analyzers  are  set  at  a  center 
frequency  of  165  kHz,  but  with  the  310A  used  as  a  low-gain 
unit.  Settings  for  the  3590A  remain  unchanged. 

310A  settings  are:  input  voltage  -  1.0  volts 

range  -  -40  dB 

The  attenuator  set  routine  now  triggers  off  channel  -Zl,  the 
low-gain  channel. 

7/8/83  -  Data-taking  ends 

7/12/83  -  Frequency  response  calibration  of  system. 
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TABLE  II 


5/13/83 

5/19/83 

5/24/83 

5/26/83 


5/27/83 

5/31/83 

6/3/83 

6/6/83 

o/ 10/83 

6/13/83 


GWEN  PROJECT 

LIST  OF  UNIQUE  CHARACTERISTICS  OF  INDIVIDUAL  TAPES 


-  Test  tape  with  one  channel  active 

-  Calibration  signal  on  rec.  #  15. 

-  Calibration  signal  recorded  on  recs.  #  39  and  40. 

-  Local  thunderstorm  observed.  A  variety  of  settings  used  to 
capture  local  noise.  404  records  recorded  before  power  fail. 

The  tape  history  is  as  follows:  Tape  started  at  13:36:00  EDT, 
just  before  rain  started  locally.  The  storm  passed  about  5  mi. 
west  at  about  13:50-13:55.  Initial  settings  were: 

input  voltage  range 

3590A  0.3  volts  -50  dB. 

310A  10  volts  -50  dB. 

No  cal.  was  taken.  After  rec.  #8,  the  310A  input  voltage  was 
changed  to  30  volts.  After  rec.  #19,  the  310A  bandwith  was 
changed  to  200  Hz,  the  3590A  bandwidth  to  100  Hz,  and  the  sample 
rate  to  1  kHz  (5  sec  sample  duration).  After  rec.  #43,  the  sample 
rate  was  changed  to  500  Hz.  After  rec.  #64,  the  310A 
input  voltage  was  set  back  to  10  volts.  There  appears  to  be  a 
~2  kHz  carrier  present  on  narrow  bandwidths.  Recording  continued 
until  rec.  #404,  when  a  power  fail  occurred  during  a  second 
thunderstorm. 

-  Data-taking  interval  set  to  5  min.  for  weekend  operation. 

650  records  recorded  before  power  fail. 

-  980  records  recorded  before  power  fail. 

-  Data  interval  set  to  36  min.  for  weekend.  200  records  recorded 
before  power  fail. 

-  The  first  part  of  this  tape  (all  data  recorded  before  02:02:40 
on  6/7/83)  is  no  good  because  of  a/d  converter  hangup.  1040 
records  recorded  before  power  fail. 

-  Data  interval  set  to  36  min.  for  weekend.  350  records  recorded 
before  power  fail. 

-  The  latter  portion  of  this  tape  (all  after  230  records)  is  no  good 
because  of  a/d  converter  hangup.  380  records  recorded  before  power 
fail . 


6/2U/83  -  120  records  recorded  before  power  fail. 


6/23/83  -  Test  tape  recorded  for  amplitude  calibration 
6/27/83  -  270  records  recordea  before  power  fail. 
6/30/83  -  Tape  stopped  manually  after  1060  records. 
7/1/83  -  Data  interval  set  to  36  min.  for  weekend. 

7/8/83  -  Log  unavailable. 


Table  III  -  Owen  project  -  tape  log 

Notes:  General  characteristics  of  data  are  annoted  in  previous  Cables. 

Uu  stop  times,  "+n“  means  that  the  tape  was  stopped  at  Lhat  time  n  days  alter  the  start  date. 
Distances  under  TV  weather  column  are  to  nearest  storm  front  from  digital  radar. 


6/10/ttJ  14:0^:50  ll;0y;5Utl  36  none  510  oti.  (t) 

1105  mi.  (WSW),  heavy 


no  report 


o  •  c 

Q.  ■£  3 

V  5 


a>  ^ 

o..  •— <  ■ 

rT  *  o 

u  0  u 

-J  a. 


3  es 
<N  -•  O 


GWEN  RECEIVER  BLOCK  DIAGRAM  AS  OF  5/19/83 


recording  equipment  setup  as  of  5/19/83. 


NEW  GWEN  RECEIVER  BLOCK  DIAGRAM 
FOR  INCREASED  DYNAMIC  RANGE  5/25/83 


figure  2  -  System  U.  Changes  as  of  5/25/63.  Two  wave  analyzers  are 
now  used  to  increase  dynamic  range. 


GWEN  RECEIVER  MODIFICATIONS  AS  OF  6/13/83 


vv 


GWEN  RECEIVER  MODIFICATIONS  AS  OF  6/28/83 


FORTH  REAL-TIME  DATA-TAKING  PROGRAM 

FLOW  CHART 


Figure  5  -  Forth  real-time  aaca-taking  program  clock  diagram. 
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MPLITUD 


GWEN  System  Response  (Test  input  to  Mpx  input) 


Figure  7  -  Measured  passband  of  noise  recording  system  from  test  input 
to  multiplexer  input.  Center  frequency  set  to  165.0  kHz. 
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APPENDIX  A 


I.  Discussion  of  Programs  Used  in  Noise  Analysis 

The  tape  which  is  written  by  the  on-site  HP-2H6C  mini-computer  consists 
of  approximately  1130  records,  each  containing  10,000  data  samples  followed  by 
a  20  word  identification  section.  The  output  voltages  from  the  two  receivers 
have  been  digitized  as  one  16-bit  word  per  data  sample  for  each  receiver. 

Since  the  two  receiver  outputs  (R1  and  R2)  have  been  multiplexed,  alternate 
data  samples  in  each  record  correspond  to  the  output  from  each  receiver. 
Successive  samples  from  the  same  channel  are  separated  in  time  by  the  sample 
rate  as  indicated  in  the  data  log.  The  A/D  converter  used  produces  an  output 
number  consisting  of  11  bits  plus  sign  (  range  -4096  to  +4096  ).  The 
right-hand  4  bits  in  each  16  bit  word  are  "garbage",  and  must  be  discarded  by 
arithmetically  shifting  each  sample  4  bits  to  the  right  (equivalent  to  a 
divide  by  16).  This  pecularity  was  not  eliminated  at  the  source  computer 
because  of  timing  constraints  caused  by  the  high  data  rate. 
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***************** 


9997  9998  9999  10000 

********************************************************** 


R1  R2  R1 


R1  R2 


R1  R2  AB  CDEFGHI JK1MNOPQRS  T 


***************** 


********************************************************** 


RECORD  N 


I.  D.  SECTION 


The  ID  section  contains  the  start  time  for  each  record  together  with  the  date, 
attenuator  setting,  time  interval  between  samples  and  the  frequency  that  was 
monitored.  Table  Al-1  (p.  65)  gives  a  detailed  summary  of  the  information 
contained  in  the  ID  record.  All  tapes  are  9-track,  800  bpi.,  with  no  pre-  or 
post-ambles  or  tape  label. 

During  the  data  processing  phase,  an  entire  data  tape  is  reaa  into  a  VAX 
11/780  disc  file  in  order  to  reduce  the  processing  time  per  tape. 


Table  Al-1. 


The  Identification  Section 


Word  Number 


Identification 


A 

Time  in  .1  second  intervals  since 
start  of  tape  (sign  bit  is  a  data  bit) 

02 

B 

Time  since  start  of  tape  in  multiples 
of  6553.6  seconds. 

C 

Year. 

04 

D 

Month. 

>05 

E 

Day  (Eastern  Daylight  Time). 

F 

Hour . 

G 

Minute . 

H 

Second. 

I 

Monitoring  frequency  in  kHz. 

10 

J 

Bandwith  in  Hz. 

>11 

K 

Location  #  (always  1). 

>12 

L 

Antenna  orientation  (always  0  degrees) 

13 

M 

Time  interval  between  groups  in 
units  of  0.1  seconds. 

14 

N 

Record  count,  starting  with  1  at 
start  of  tape. 

0 

#  of  records  in  group. 

>16 

P 

Attenuation  setting  used  (1,2,4  or  8). 

A.  PROGRAM  NOISE 


This  program  reads  each  record  from  Che  VAX  data  file  chrough  Che  code 
combinacion: 

CHARACTER  BUFFER*32768 
BYTE  JDATAC 32768) 

EQUIVALENCE  (JDATA, BUFFER) 

**** 

**** 

OPEN  (UNIT=1 , CARRIAGEC0NTKOL= ' LIST' , STATUS® ' UNKNOWN 1 , RECORDS IZE*20040) 


READ (1,5, END* 1004)  BUFFER 
FORMAT (A) 


CLOSE(UNIT=l) 

An  encire  record,  consisCing  of  20040  byCes,  is  read  inCo  an  array  called 
BUFFER.  The  HP-2116C  wriCes  a  2-byCe  daca  word  onCo  Cape  wich  Che  lease 
significanc  byce  firsc.  Since  Che  VAX  11/780  cape  I/O  sofeware  expects  Co  see 
Che  mosC  significanc  byCe  firsC,  che  byCe  order  in  each  of  Che  10020,  2-byCe 
words  musC  be  reversed.  This  operacion  is  accomplished  by  a  call  Co  Che 
subrouCine  DECODE.  Here,  Che  byce  order  is  swapped  using  a  MACRO  routine 
called  BYSWAP.  The  resulting  2-byte  word  is  then  passed  back  to  the  main 
program  through  the  INTEGER*2  array  JDATA.  In  addition,  the  contents  of  the 
ID  record  are  stored  in  the  array  IDREC  ana  the  record  date  and  start  time  is 
converted  into  'year,  month,  day,  hour,  minute,  second'  format. 

The  contents  of  the  first  100  daca  samples  are  tested  against  the 
corresponding  data  words  in  the  previous  record  to  insure  that  record 
duplication  has  not  occurred.  This  would  be  the  case  if  the  multiplexer 
malfunctioned.  The  current  record  would  then  be  ignored  in  subsequent  data 
analysis.  For  the  first  record  of  each  tape,  the  corresponding  power  level 
from  each  receiver  output  is  determined  by  squaring  the  digital  value  of  each 
data  sample  in  Che  record  and  multiplying  by  a  calibration  constant.  This 
background  power,  in  units  of  10”^W/m^,  is  then  plotted  as  a  function  of 
time  for  Che  first  record. 

The  statistical  characteristics  of  a  particular  record  are  determined  in  a 
subroutine  called  STAT.  The  Digitally-Controlled  Attenuator  has  4  states 
(1,  2,  4,  8  in  ID  word  ' P 1 )  and  the  ID  record  identifies  which  of  the  4  is 
operative  at  the  time  a  record  was  taken.  The  attenuation  is  compensated  for 
in  Che  statistical  analysis  by  dividing  the  digitized  voltages  in  each  record 
by  the  factors  corresponding  to  the  appropriate  attenuation  step,  namely  1.0, 
.0631,  .00668,  .000708.  The  resulting  values  are  then  squared  ana  multiplied 
by  the  calibration  constant  appropriate  to  each  receiver,  measured  at  the 
attenuation  setting  of  1.0  to  determine  the  instantaneous  power.  The 
quantities  that  are  determined  for  each  record  and  for  each  receiver  output 
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For  che  data  taken  after  June  15,  1983,  only  Pm  t  the  highest  noise  spike 
found  in  a  record,  was  determined.  The  quantity  is  a  measure  of  the 
impulsiveness  of  the  noise  in  a  given  record.  The  quantities  (,  P>  ,  VD  and 
Pm  are  thereby  determined  for  each  record  and  plotted  as  a  function  of  local 
civil  time. 

B.  Program  APD 

The  Amplitude  Probability  Distribution  (APD)  is  determined  by  reading  a 
data  record,  compensating  for  any  attenuation  changes,  and  computing  the 
percentage  of  time  in  a  single  record  that  the  background  noise  exceeded  a 
prescribed  threshold.  The  average  background  power  level'CP^of  each  receiver 
is  computed  and  a  total  of  100  threshold  levels  are  determined  by  simple 
multiples  of  <P>.  The  number  of  data  samples  exceeding  each  threshold  is  then 
found.  The  result  is  then  plotted,  for  a  single  record,  on  a  Log-Log  graph. 
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c. 


Program  TPD 


Once  a  particular  record  has  been  selectee  from  the  data  file,  a  Time 
Probability  Distribution  (TPD)  is  computed  as  follows:  The  detected 
background  power  P  for  each  data  sample  is  corrected  for  the  attenuation  as  in 
previous  programs.  A  threshold  power  level,  T0>  is  selected  and  a  function 
P  is  defined  such  that, 


Pi  -  1 

for 

Pi  >  T0 

*  0 

for 

Pi<T0 

where  i  is  the  index  of  the  data  sample  in  each  record.  Given  this  function 
pi,  the  autocorrelation  function  A(j)  is  computed 

5000 

E  ?i  Pi+j 

i  =  l 

A(j)  =  _ 

5000 

E  ^  n 

i*l 


where  j  corresponds  to  the  offset  period  t  in  the  integral  definition  of  the 
autocorrelation, 


A(  t  ) 


00 

P(t)*P(t-T  )dt 


This  function  is  a  measure  of  the  fraction  of  noise  pulses  above  the  threshold 
power  T0  that  have  a  separation  in  time  of  j  x  .16  milliseconds.  The  TPD  is 
related  to  the  autocorrelation  function  A(j)  by 


V 


T(  >  k  )  -  j  >  k 


5000 

E 

j-i 


A(j) 


A(j) 


where  A(j)  is  the  total  percentage  of  pulses  above  the  threshold 
and  che  quantity  that  is  actually  plotted  on  the  Log-Log  vs  Log  graph  is 

1 

-Logi0(Logi.o( - ))  . 
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3  10U00  CONSTANT  LBUF  LBUF  i  LI  I  /  CONS 'I  ANT  NCO 

4  LUt.ir  20  f  OOIM  OI'.MJF  0  0 1  Ul  I.BUF  +  CONSTANT  ID  BLOCK 

5  2  L  0  LOAD  (  ’? -TRACK  >  271  LOAD  (  MINI  VENTER  DRIVER  ) 

6  0  f).UUF  LAD  !  LBUF  20  NLQC  ! 

7  275  LOAD  (  ID  WORDS  ) 

8  272  LOAD  3  FLO  ! 

?  L 7 R RESET 

10  ;S 
1  1 
12 

13 

14 

15 


BLOCK  271  (417EO 

0  (  MINIVER  TER  DRIVER  FOR  GWLN  WITH  AUTOCYCLE  ) 

1  BASE  9  OCTAL  17  CONSTANT  MX  22  CONSTANT  OBSLOT 

2  0  VARIABLE  MAD  0  VARIABLE  NSP 

3  CODE  SAMP  S  )  0  LD,  0  TC,  NSP  0  ST,  SI  )  0  LD,  MAD  0  ST, 

4  20000  *  0  LD,  MX  OTA,  MX  STC ,  ,C  ELGIN  MX  SFS,  END 

5  ELGIN  0  CLF,  OBSLOT  STC,  ,C  BEGIN  OBSLOT  SFS,  END 

6  OBSLOT  CLC ,  ,C 

7  1  00017  I  0  L.D ,  7  OTA,  MAD  0  L.D ,  1  0  0000  *  IOR,  3  CLC, 

S  3  OTA,  #CH  MINUS  *  0  LD ,  3  STC,  3  OTA,  40000  *  0  LD,  MX  OTA. 

9  MX  STC,  ,C  7  STC,  ,C  0  STF ,  MAD  0  LD ,  *CH  *  0  AD,  MAD  0  ST, 

10  NSP  ISZ,  END  2P0P  , 

11  7  ICQDE  UP LX  7  CLF,  JMP ,  BASE  ! 

12  ;S 

13 

14 

15 

BLOCK  272  (420B) 

0  (  GWEN  COVER  PROGRAM  ) 

1  10  VARIABLE  NRLC 

2  !  FILLARR  0  OBUF  NSQ  SAMP  i 

3  276  LOAD  273  LOAD 

4  :  DOIT  NR  EC  (3  0  DO  FILLARR  "('WRITE  LOOP  ; 

5  i  WAIT  400  0  DO  LOOP  ;  (  ABOUT  20  MS,  ) 

6  :  OGLOIS  TOD  8  IDBLOCK  !  TOD  1+  5  IDBLOCK  1+  !  IDBLOCK  14  +  @ 

7  0  DO  FILLARR  I  WRITE  WAIT  LOOP  ; 

8  :  PR PAR AM  CR  20  0  DO  IDBLOCK  I  +  S  S ,  LOOP  ; 

9  :  I  NT  TEST  BEGIN  TOD  @  IDBLOCK  9  -  DUP  0<  IF  32768  +  1 IIEN 

10  IDBLOCK  12  +  B  1  -  >  IF  1  ELSE  GTLST  0  THEN  END  j 

11  :  DO  I  APE  IDSET  REW  GAP  1  0000  0  DO  INTTEST 

12  DOLOTS  PRPARAM  FINDSTATS  EOT?  IF  WI'M  RWO  ABORT  THEN  luJI’  ; 

13  s& 

14 

15 
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t 


LMJ'II'U  .'75  \  un 

o  (  .'i/O  a <  o- : . i • ,  n  1 1 k  >  .  .  ■  ns  > 

1  :  a/ d dump  c  o  o  no  t  2  x  s ...  jr  ;•  i6  /  .  i  2  x  t  ♦  u t:<ur  a 

2  16/  .  cR  I  ‘I UP  ; 

.5  1  000  VARI..‘'LE  ILSTAttP 

4  :  3CAUARRA  (  0  NGQ  0  DO  l  2  x  unUf  ft  ADS  I  EG  I  AMP  ft  >  II  0-  LEAVL 

5  l  lll  .N  LOOP  ; 

6  :  WAIIl  URBUPST  ELGIN  F  Tl  I  ,*kR  ■  LANARK  AY  END  ; 

7  :  ALARM  REW  1 1)  GET  1  0000  0  DO  AA  L  I  i  UK  EUR  ST  TWRIFE  EOT?  IF 

3  WIM  R WO  ABORT  THEN  LOOP  ; 

?  0  OAR  I  ABLE  MAX  0  VAR  I  ABLE  f'IN  0,  REAL  AVG 

10  ;  FINDS  TATS  2  0  DO  0  MAX  !  0  MIN  !  0.  AVG  F!  500  0  DU 

11  I  2  x  J  +  O OUF  ft  16  /  DUP  MAX  ft  >  IF  OUP  MAX  ! 

12  ELSE  DUP  MTN  9  <  IF  DUP  MIN  •  THEN  THEN  SFLOA1  FSG  AVG  F+ ! 

13  LOOP  CR  I  .  AVG  Fft  NSG  SI  LOAF  F/  FSOkT  F,  MAX  ft  S.  MIN  ft 

14  S,  LOOP  ; 

15  > 

BLOCK  274  (42210 

0  (  MORE  GWEN  GOODIES  ) 

1  :  PWRRLSET  SCL  CLR  3  TBGSTART  1  ATTSET  }, 

2  ;  S 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

BLOC 
0 
1 
2 

3 

4 

5 

6  .  11  BW?"  SASK  IDBLOCK  9  +  !  LOCATION?"  SACK  IDBLOCK  10  +  ! 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 


40 


ft 


275  (  423 B ) 

ID  BLOCK  SETUP  FOR  GWEN  DATA  ) 

TSET  . "  YR?"  SASK  IDBLOCK  2  +  !' 

'  MONTH?  *  SASK  IDBLOCK  3  i-  !  .  "  DAY?"  SASK  IDBLOCK  4  + 

1  ||R?"  SASK  IDBLOCK  5  +  !  .  "  MIN?"  SASK  IDBLOCK  6  +  ! 

•  SEC?"  SASK  IDBLOCK  7  <  !  0  TOD  !  0  TOD  1+  T  j 

LOGSET  FRLQ? "  SASK  IDBLOCK  8  +  ! 

1  BW?"  SASK  IDBLOCK  9  +  !  . "  LOCATION?"  SASK  IDBLOCK  10 

1  ANT.  ORIENTATION?"  SASK  IDBLOCK  11  +  ! 

1  SAMPLING  INTERVAL?"  SASK  IDBLOCK  12  +  \ 

•  *  OF  RECORDS?"  SASK  IDBLOCK  1 4  <■  !  •, 

IDSET  TSET  LOGSET  16  IDBLOCK  +  4  AZERG  ; 


-  .*•  A 


>VV 


-IJ. 'M  ii-'t,  1  •  1 1;  > 

0  (  }>  i  C !  I  ,  .1  L.'t  i  i a  I  1 .1 1 1  I  L.0  A  I  I  •  : '  I  '.'I  K  ) 

1  . (I  b  '  i..i  r  • '  •  i  .'if!  I  >||  ICH  1  .-.c  1  C  ,M  IVAI  ml  10  0  t'Dlri  «■_.  r 

2  u.ioori  i;( •i-.'oi ant  dm, a, mc  o  vaui.pil  dovnir 

3  CUD::  AffSEI  S  )  o  LD,  ill  ICH  IMA.  POP  . 

4  :  APR  l  t:;R  .  ••  ATT.  SETTING  C  | ;  AH' 11:  D  10  M  AlfVAliJL  ?  ; 

5  :  T  A  (  0 ET  FOB  @  DOWNCT  !  ; 

6  ;  AFI  CLT  DUP  AT  I  GI..T  A  I  TWO  Ul  *  ;  1  A I  I  SET 

7  :  A  HUP  A  RVALUE  &’  DUI“  0  <  IF  2  *  AlfGLT  APRT 

8  TAT  GET  LLSL  DROP  THEN  ; 

9  :  A  IT  DOWN  AH  VALUE  8  DUP  1  >  IF  2  /  A  IT  SET  APRT 

10  1  A  T  S E T  £ L. S E  D POP  THEN  ;  ,(  t 

11  1024  CONSTANT  LJPTEST  200  CONSTANT  DOWN  TEST 

12  :  UP  I  LG  0  5  0  0  DO  I  2  x  DETEST  6?  ADS  UP  TEST  >  IT  0=  LEAVE  THEN 

13  LOOP  ; 

14  --> 

15 

BLOCK  277  <42SB> 

0  (  hUPE  DIGITALLY  CONTROLLED  ATT.  DRIVER  > 

1  :  DOUNFLG  1  50  0  DO  I  2  »  DETEST  3  AES  DOWNIEST  >  IF  0= 

2  LEAVE  THEN  LODI*  5 

3  i  DO  TEST  0  DETEST  50  SAMP  j 

4  :  DOTEST  DOTLST  100  0  DO  I  DETEST  @16/1  DETEST  !  LOOP  i 

5  :  LAZY  TOD  @  DOWNCT  @  -  DUP  0<  IF  32768  +  THEN  DWNINT  > 

6  IF  1  ILSE  0  THEN  ; 

7  :  Li  TEST  DOTEST  UPFLG  IF  AT  TUP  ELSE  DOWNFLG  IF  LAZY  IF  ATTDOWN 

8  THEN  THEN  THEN  ATTVALUE  @  IDELOCK  15  +  !  t 

9  •>  S 
10 

1 1 
12 

13 

14 

15 

BLOCK  278  (426B) 

0  (  MIN1VERTER  DRIVER  FOR  GWEN  ) 

1  EASE  6  OCTAL  2  CONSTANT  iFCH  17  CONSTANT  HX  22  CONSTANT  DESLOT 

2  0  VARIABLE  HAD  0  VARIABLE  NSP 

3  CODE  SAMP  S  )  0  LD,  0  1C,  NSP  0  ST,  SI  >  0  LD,  MAD  0  ST. 

4  ELGIN  0  CLF  >  OBSLOT  STC,  ,C  BEGIN  OBSLOT  SFS,  END 

5  OBSLOT  CLC ,  .C  20000  $  0  LD,  MX  Ol'A,  MX  STC,  .C  HERE  MX  SFS. 

6  IMP.  10  0  017  *  0  LD,  7  OTA.  MAD  0  LD.  1  0  0  0  0  0  *  I  OR,  3  CLC. 

7  3  OTA.  tfCH  MINUS  T  0  LD,  3  STC,  3  OTA,  40000  *  0  LD .  MX  OTA, 

8  MX  STC,  ,C  7  STC,  ,C  0  STF,  MAD  0  LD ,  D=CH  1  0  AD,  MAD  0  ST, 

9  NSP  ISZ,  END  2P0P  , 

10  7  ICODE  UPEX  7  CLF,  JMP,  BASE  ! 

1 1  ;S 

12 

13 

14 

15 


ft  .. 
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.*■■  y-VV-V.VV- 


i  ;  ,  -  .  it  > 

0 

1 

■  > 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

BLOCK  280  (  430  E< ) 

0  (  DIGIDATA  9-TRACK  ROUTINES  )  BASE  6  OCTAL 

1  23  CONSTANT  MID  MTD  1+  CONSTANT  MTC 

2  23  VARIABLE  RCD  31  VARIABLE  WCD  0- VARIABLE  BOB 

3  CODE  STATUS  MTC  LIA,  PUSH  , 

4  :  BUSY?  BEGIN  STATUS  1410  AND  0=  END  j 

5  ;  PE?  STATUS  2  AND  2  /  ; 

6  :  REJ?  STATUS  10  AND  10  /  j  ;  EOT?  STATUS  40  AND  40  /  ; 

7  :  DOT?  STATUS  100  AND  100  /  j  :  EOF?  STATUS  200  AND  200  /  j 
Q  :  RING?  STATUS  4  AND  4  /  1  XOR  j 

9  :  LOCAL?  STATUS  1  AND  ;  :  TRANSB?  STATUS  1000  AND 

I  0  1  0  00  /  ;  •.  RLW?  STATUS  2000  AND  20  00  /  ; 

II  CODE  CMMD  S  >  0  LD,  MTC  OTA,  MTC  SlC,  ,C  POP  , 

12  CODE  MTB  MTC  CLC,  ,C  MTD  CLC ,  ,C  1510  11  0  LD ,  MTC  OTA, 

13  NEXT  , 

14  — > 

1 5 


BLOCK  281  <  43 1  EO 

0  <  MORE  9-TRACK  ) 

1  :  GAP  15  CMMD  BUSY?  ;  :  lrSR  3  CMMD  BUSY?  ;  :  GSR  41  CMMI) 

2  BUSY?  ;  :  RWO  105  CMMD  ;  :  REW  101  CMMD  BUSY?  ;  s  CL R  110  CMMD 

3  :  SEL  1400  CMMD  ;  :  WF M  211  CMMD  BUSY?  j 

4  !  FSF  203  CMMD  BUSY?  ;  :  JBSF  241  CMMD  BUSY?  i  :  GFM  215  CMMD 

5  BUSY?  j  SEL  CLR 

6  BUF  9  VARIABLE  SAD  1001  VARIABLE  NLQC 

7  CODE  TREAD  MTD  *  0  LD,  6  OTA,  2  CLC,  1  0  0  0  0  0  H  0  L.D ,  SAD  0  AD, 

8  2  OTA,  2  STC ,  NLQC  0  LD ,  0  TC,  2  OTA,  RCD  0  ID,  BOB  IOR ,  MTC 

9  OTA,  MTC  STC,  ,C  MTD  SIC,  ,C  6  SIC,  ,C  NEXT  , 

10  CODE  TWRITE  MID  ♦  0  LD,  6  OTA,  2  CLC,  SAD  0  LD ,  2  OTA, 

11  2  STC,  NLOC  0  L.D  ,  0  TC,  2  OTA,  WCD  0  LD,  BOB  IOR,  MTC  OTA,  MIC 

12  STC,  ,C  MTD  STC,  ,C  6  SIC,  ,C  NEXT  , 

13  6  ! CODE  TP  I  6  CLC,  MTD  CLC,  ,C  MTC  CLC ,  JMP  , 

14  1  +BLOCK  CONTINUED 

15 
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ll 

1 

'■J 

3 

4 

5 

6 

7 

8 
9 

1  0 
1 1 
12 

13 

14 

15 


:  i  1 1  I  i  .  N  1  I  I-  1  I  ’  M  .  j  ,  ,  | 

;  I  ••Kilt  UK.,.’  i  < 1  v 
:  POO  I'l'O  !)  IKL.C  ! 
ItHUd  NI  OC  !  Il  l:  AD 
10 1)1  NLOC  1  1URI.  i  E 


t 

v  Oil  i.  i.  i  1 1  * i  ;  '  • 

f)|  Cl.'iAL  U‘CI  i  ILK  I  > 

1  l-'l.AD  BUSlM  IkI.D  ; 

KM.  W  k'LW  0  i  l-M  C  i  ; 
t:  IK  LA  I)  BUI  It?  MAI)  ! 

B  I  KIR  J.  IE  BU:  (t?  SAD  ! 
iBl..K  PREV  (>?  !  ; 

BIN  0  BOB  !  :  BCD  id  mi  .in  ftflrt  i 

RE?  PE?  IP  7  Erl  IT  . 

s: 

BNO  pt:?  Pku::v  (•?  ?  status  o.  \ 

MSP  DUP  0<  IF  MINUS  0  DO  BCR  LOOP 
i'iCr  DUP  0<  II"  MINUS  0  DO  BSF  LOOP 
BASE  !  ;  S 


i  1 


1  0  0  0  0  0  BOB 
>1  -  1 1  SEN 


BSF  BSR 

ELSE  ( 
ELSE  C 


:%■  re  *  ! 


DO 

DU 


F  SR 
FSF 


LOOP 

LOOP 


BLOCK  203  <43310 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 


BLOCK  (434EO 

0 
1 
2 

3 

4 

5 

6 
7 
Q 
V 

10 
1 1 
12 

13 

14 

15 
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THEN 

THEN 


■LUCK,  '  CL,f:) 

o  <  i  i  mi.  call  gt  ui  pmor  i ■  <  *: ri  i,jls  ) 

1  CAGE  U  i)|jp  OCTAL  ion  COi;:.,  TANI  I  LI  I 

2  I"  i  ND  !  C  (.  IDE  0  ~  1 1  1  R  U  f  \.  Nil.  :■  K  U  P  I  l  I  1.  N  D 

3  CODE  r  DCS  I  UP  (CM  CLC ,  NeXT  . 

A  I BGSTOP 

5  CO  IDE.  1  OUST  ART  TCI!  CLC,  S  )  0  LD ,  TCI!  I>UP  01  A,  BTC,  ,C  POP 

6  0  VARIABLE  TOD  ,  0  VARIABLE  ELI 

7  TCH  !  CODE  TOC  INI  TCH  CLF  ,  TOD  IBZ ,  DLIP  DUP  JMP  ,  TOD  It  IBZ 

8  JT1P ,  JT1P  , 

9  :  SWB  TOD  &  ELT  !  ;  ,i  SWP  TOD  8  ELI  8  -  .  ; 

10  BASE  ! 

11  ;S 

12  TOO  STARTED  BY  N  TBGS1ART 

13  INTERVAL  IS  1  SEC  /  <10  A  [4NI  > 

14  I.E.,  A  IS  1  SEC,  1  IS  1  MILL I SEC ,  ETC. 

15 

<LQCI<  286  <436B> 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


LOCK  287  <437B> 


APPENDIX  B  -  PROGRAM  LISTINGS 

7i\  is  i*  4  listing  of  the  actual  FORTRAN  program  ;od«s  used 
.'■nylementing  the  tap#  reading  and  data  analysis  operations 

PROGRAM  ngise 

THIS  PROGRAM  *AS  WRITTEN  3V  DR  STEN  GDENWAlD  FOR 
THE  GWEN  PROGRAM  AT  NRL  ON  MAY  R,  1=83 

FUNCTION:  TO  READ  A  9-TRACK  TAPE  AND  EXAMINE  “HE 
NOISE  STATISTICS  OF  THE  VLF  RADIO  BACKGROUND 

THIS  PROCRAM  IS  SPECIFICALLY  DESIGNED  TO  HANDLE  THE 
TAPES  RECORDED  AFTER  MAY  31.  1983  WHICH  CONTAINED 

DATA  GROUPS  CONSISTING  OF  10  CONSECUTIVE  DATA  RECORDS 
THE  PROGRAM  CALCULATES  AVERAGED  QUANTITIES  FOR  EACH  GROUP 
AND  CAN  BE  MODIFIED  T0  HANDLE  THE  DATA  BETWEEN  MAY  20  AND 
MAY  31 


COMMON/HEADER/  IDREC<30>.  1ST  <  4  J , STSC, SPSC 
COMMON/DATA/  I DATA! 2, 10000 > , NCH 
COMMON/AUTO/  TPD < 2. 3. 5000 ) 

COMMON/TIME/  IH.  IM, SEC.  ATIM<  5000) 

CHAR AC TER *30  FILENAM 
CHARACTER  BUFFER *32768 
BYTE  JDATA <32768),  TEST ( 100) 

EQUIVALENCE  (JDATA.  BUFFER) 

DIMENSION  FMAXC2, 2000) 

DIMENSION  TPDI (3). X<5000>. PERCENT ( 2. 5000). SIGMA(2. 2000) 
DIMENSION  AVE ( 2. 2000).  FTIM<5000), FDAT A ( 2.  5000/.  TIM(2000> 

LOGICAL  ERRFLG 

OPEN  THE  TAPE  FILE  ON  THE  SCRATCH  DISC  CONTAINING 
THE  DATA  FROM  THE  GWEN  OBSERVING  RUN 

TYPE*.  'INPUT  TAPE  FILE  NAME' 

READ <5. 333)  FILENAM 
FILENAM-'SIAl.  C0DENWALD3B  DAT  • 

3  FORMAT < 5X.  A30) 

OPEN (UNIT*!,  CARR  I AGECONTROL=  'LIST  ' . 

1  STATUS- ' UNKNOWN ' ,  RECORDSI ZE-20040) 

I DATA  CONTAINS  THE  INTEGER  VALUES  FOR  UP  TO 
10.000  DATA  SAMPLES  FOR  TWO  CHANNELS  EACH 
HAVING  N1  AND  N2  SAMPLES  RESPECTIVELY. 

I TPD  IS  THE  VECTOR  CONTAINING  THE  THREE 
THRESHOLD  SENSITIVITY  LEVELS  FOR  THE  TPD 


PROGRAM  BEGINS 

DEFINE  NUMBER  OF  WORDS  IN  DATA  RECORD 

ID  RECORD 

NDA  -  10000 
NID  -  20 
NTOT-NDA+NID 
NBY  -  2*NT0T 
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NUMBER  of  tape  records 


MREC  »  150 

nskip  =  o 

DO  38  1  =  1. 30 
IDREC*  I  ■  =0 
ib  CONTINUE 

calioration  constant  for  gain  1  to  convert 

»ro«  data  units  (volts)  squared  to  watts/meter  square 

P c a  1 =8  62e-l? 


initialize  plot  scales 


<S1  =  002 

*S1  =  1  OE-6 

uOG-LOG  plot 

<S2  =  3.  33 

YS2  =  1.33 

AVERAGES  AND  SIGMA 

XS3  =  005 

VS4  =  2.  OE-12 

YS5  =  005 

YS3  -  002 

TYPE*,  'TAPE  NUMBER  ' 

READ*,  IDREC  <  30  > 

TYPE*,  'NUMBER  OF  RECORDS  TO  SKIP' 

READ*, NSKIP 

type*,  'NUMBER  OF  RECORDS  TO  PROCESS' 
read*, mrec 

TYPE*,  'SKIP  =  ', NSKIP,  '  PROCESS  «  ' , MREC 
IDREC  <30)=1 


INITIALIZE  PLOT  UTILITY 
CALL  PLOTST  (  5,  'IN',  0) 


C  INITIALIZE  DATA  ARRAYS 

DO  35  1*1,  2 
DO  40  J-l.  100 
PERCENT ( I,  J)=0 
DO  41  K-l,  3 
TPD (  I,  K,  J )  =0 
Al  CONTINUE 

AO  CONTINUE 

35  CONTINUE 

DO  7 fa  1*1,  100 
TEST ( I  )  «  0 
CONTINUE 
C 

C  READ  TAPE  RECORD  OF  LENGTH  NBY  BYTES 

ISTT-1 
IPAR  *  0 
NCH  »  2 
IPLOT  *  1 
NNREC -NSKIP 
NREC  *  0 

IF1NSKIP  EGO)  GO  TO  100 
DO  201  11=1, NSKIP 
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READd.5)  BUFFER 
FORMAT (A) 

.'.'1  CONTINUE 

MXREC*MR£C+NSKIP 
WRITEC6,  555) 

;fr  fqamaT (  '  TAPE  RECORDS  SKIPPED') 

FORMAT ( 1HO. 6X .  'TIME'.6X.  '  RECORD',  8«.  CHANNEL'.  6X, 
l  AVERAGE '.6X,  '  v-SUB-D/ 20DB  ' .  10X,  'PEAK/M  SO'///) 

IFIRST  =  1 
jgain  =  0 

10  NREC*NREC*1 

IF < NREC  GT  MRECi  GO  TO  1005 
READ (1.5.  END= 1 004 )  BUFFER 
WRITE ( 6, 888 )  < JDATA< I ).  1*1.  10) 

3  586  FORMAT (1015) 

NOW  TEST  THE  FIRST  100  BYTES  TO  MAKE  CERTAIN 
0  that  THE  CURRENT  RECORD  AND  THE  PREVIOUS 

I  RECORDS  ARE  NOT  IDENTICAL.  IF  THE  MULTIPLEXER 

HANGS  UP  THIS  WILL  BE  THE  CASE 

IF< IFIRST  EG. 0)  GO  TO  S89 
C  THIS  IS  THE  FIRST  RECORD  FROM  THE  TAPE 

DO  87  IM*1, 100 
TEST  < I >*JDATA< I ) 

37  CONTINUE 

IFIRST  *  0 
GO  TO  887 

389  DO  86  1*1. 100 

IF  (TEST  < I ) .  NE  UDATA(D)  CO  TO  887 
36  CONTINUE 

IF  ALL  ELEMENTS  NOT  IDENTICAL  THIS  IS  GOOD  DATA 
WRITE (6. 886)  NREC 

386  FORMAT (  '  DETECTED  REPEATED  RECORD  '.14.  '  RECORD  SKIPPED') 
NREC*NREC-1 
GO  TO  100 

567  CALL  DECODE! JDATA, NBY. NREC. I STT) 


C  see  if  the  gain  setting  has  changed 

IF ( JGA IN.  EG.  IDREC ( 16) )  GO  TO  937 
WRITE! 6. 938)  IDREC(16) 

538  FORMAT  < '  THE  GAIN  IS  NOW  SET  TO  15) 

JGA IN* IDREC ( 16 ) 


COMPUTE  AVERAGE  VALUE  AND  SIGMA  OF 
DATA  IN  RECORD 

PLOT  THE  FIRST  RECORD  IN  CHANNEL  l  AND  2 
IF  THIS  IS  THE  FIRST  RECORD  PROCESSED 
EXPRESS  X  AXIS  AS  TIME  IN  MILLI  SECONDS 


°37  IFdPLOT  NE.  1)  GO  TO  50 

DELT-1000  /6200 
SUM 1*0 
SUM2*0 
KK*1 

DO  32  1*1.5000 
FTIMd  >*I*DELT 

f data ( 1 . i ) * ( 1  0* i da ta < 1 . i ) /16  )**2 
fdataf  2. l >»< 1  0*1 DATA <2.  I >/ 16  >**2 
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SUM1*SUM1«-FDATAi  2,  I  '• 
SUM2«SUM2+FDATA!2.  I ) 

CONTINUE 
AVEKK* SUM1/ 5000 
AVELL-SLM2/5000 
XST-10.  /<  5000  *D£LT. 

VS1-  2 /AVEKK 
VSO  *  2/AVELl 

WRITE! bi 89)  ( I DREC  <  > . *KK=1.  20 > 

FORMAT!'  TAPE  ID  20110/) 

IFtNREC  EQ.  1)  WRITE(6.  Ill) 


determine  the  g  axis  limit  IN  UNITS  OF 
10! -1 5 )  WATTS/METER  SQUARE/HZ 


one  channel  is  at  a  gain  of  1  jih  1 1  e  the 

second  channel  is  30dh  lauf  m  jam  <factoT  of  33) 


SCALE* VS1 

IFIFDATA! 1.  10). LT  FDATA12,  10))  SCALE-YSO 
B88-8.  *pca 1*1  Oe 1 5* ! 10** ! 30.  /20  > ) /SCALE 
AAA-8.  *pcal*l  0* 1  5/SCALE 

cl  *  aaa 

c2  *  bbb 

if(fdata<2.  1).  gt.  fdatall.  1))  cl  *  bbb 
i  f ! f da ta!2,  1 ) .  g t  fdata(l,l>)  c2  *  aaa 
Cl-Cl/3100. 

C2-C2/3100. 

CALL  PLOT!.  5.  2  .  -3) 

CALL  PLOTT <FTIM.  FDATA,  3.  XST, SCALE.  5000.  1, cl. c2) 

IPLOT-O 

DO  150  J-1.2 

CALL  STATIJ,  SIC,  AU,  AAV.  10000. NREC. FAX ) 

SIGMA! J. NREC) =AAV 
AVE(J.  NREO-AV 
FMAXU.  NREC) -FAX 

FDATA! 2, NREC)*! 1  0*IDATA(2, 1000) ) **2 
CONTINUE 
GO  TO  100 
WR ITE !  6,  2 ) 

FORMAT!'  ERROR  IN  INITIALIZING  TAPE') 

GO  TO  1000 
WRITE! 6, 3) 

FORMAT!'  END  OF  TAPE  ASSUMED') 

GO  TO  1003 
WRITE (6, 1006)  NREC 

FORMAT!'  END  OF  DISH  OR  TAPE  FILE  DETECTED 
AFTER  RECORD  NUMBER  .15) 

NREC  *  NREC— 1 


NOW  COMPUTE  AVERAGES  FOR  EACH  GROUP  OF  TEN 
RECORDS 


SMTi-O. 

SMT2-0. 

DO  412  K«  1 . 2 
J-l 

AT-ATIMI1 ) 
SUM 1-0. 
SUM2-0 
SMX*  0 


y.'V.'.i 
*•  V.' 

.V  .1 


•  V  V  1 
v-.-v 
«.  *.  1 
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DO  41 1  1*1,  NREC 

IF  ( ABS  <  AT  I  li  <  I  )  -  AT )  ST  01)  GO  TO  413 
sumi«sumi*ave;k,  i  > 

SUM2«SUM2+SIGMA;K,  I  ) 

IF(FMAX(K,  I)  GT  SMAX ) SMX  =  FMAX ( K.  I ) 

IF  <  I  EG  NR EC i  GO  TO  413 
GO  TO  411 
AVECK, J)»SUM1/10 
3 I GMA 1 *SUM2/ 10 

3IGMA<  A.  j  )  »20*AL0G  1 0  (  AVE  <  K,  .j )  /SIGIiAl  , 

FMAX ( K ,  J ) *SMX 
FDATA(  1, J ) *AT I M ( 1-1 ) 

FDATAC2. J)=FDATA(2,  I - 1 > 

AT«ATIM< I ) 

SUM  1  -A  V£  ( K .  I  ) 

SUM2=S I GMA  <  K,  I  ) 

SMX  =*  CMAX (K.  I  ) 

IF  ( K.  EG  1)  SMT 1  =  SMT1--AVE;  t.  J) 

IF(K. EQ. 2)  SMT2  =  SMT2  +AVE ( 2.  J ) 

J*U+1 

CONTINUE 

J«U-1 

SMT 1  =  SMT 1 / J 
SMT2=SMT2/U 
AVE1K.  J )=SUM1 / 10 

SIGMA ( X, J) “20  *ALOG 10 ( AVE1 K, J  > /SIGMA1 ) 

FMAX (K.  J)*SMX 
FDATA< 1, J)=AT 
CONTINUE 

WRITE  <7,  434)  IDREC  (  5  )  j  IDRECC4).  IDRECO),  1ST  (  1  ) .  IST(2) 
FORMAT('  '.514) 

DO  415  I-l.  J 

TYPE*.  I. FMAX ( 1.  I ), FMAX <2,  I ) . S I GMA ( 1,1), SIGMA (2,  I ) 
WRITE; 7, 435)  I. FOATA< 1,  I ) . AVE ( 1,1), AVE <2,  I ) 

FORMAT; '  ',14. 3F1S  2) 

atim; I )=fdata; i,  i > 

CONTINUE 
XS3  -  10.  /30 

COMPUTE  Y  AXIS  SCALES  AND  LIMITS 
YS1  -  1.  /SMT 1 

if;ave;i,  10)  lt  ave<2, io>>  ysi  =  1. /SMT2 

B8B=8  *3.  08E-9/YSI 

AAA«8  *2.  02E-4/YS1 

TYPE*.  YSI.  BBB. AAA 

Cl  »  BBB 

C2  »  AAA 

IF<AVE( 1.  10)  LT.  AVE12,  10) >  C2  =  BBB 
IF(AVE< 1.  10)  LT.  AVE12,  10) )  Cl  =  AAA 
CALL  PLOT; 11  .  0.  ,  -3) 

CALL  PL0TT1ATIM,  AVE,  5.  XS3, YSI, J, J,  cl, c2> 

PLOT  THE  V  SUB  0  FOR  EACH  RECORD 

YS5  -  4  0/100 

CALL  PLOT; 11.0,0.,  -3) 

CALL  PLOTTfATIM, SIGMA,  4,  XS3,  VS5, J,  J,  100.  .  100  ) 

PLOT  the  PEAK  POWER  IN  EACH  RECORD  '  MEAN  POWER 


CALL  PLOT ( 11  .0  .  -3) 

SC ALE-4  /l  OEll 
AAA” 06 2 

DBB-862  *10*+ ( 30  /20  ) 

C2-AAA 

Cl-BSB 

IF ( FMAX (1.1)  GT  FMAX (2. 1 ) )  Cl-AAA 
IF ( FMAX (1.1)  GT  FMAX  <2,  1 ) )  C2-BBB 
CALL  PLOTT ( ATIM.  FMAX,  6,  XS3, SCALE, J,  J,  cl, c2> 

now  plot  tape  id  parameters  on  plot  header 

I DREC  <  29 ) “NREC 
IDREC ( 28)  -  IPAR 
I  DREC ( 27)  »  1 
CALL  PLOT (9  5.  0,-3) 

CALL  HEADER (NID) 

CALL  PLOTND 
CLOSE ( UNI T-l i 
STOP 
END 
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PROGRAM  PEAKS 


THIS  PROGRAM  WAS  WRITTEN  BY  DR  STEN  ODENWALD  -OR 
THE  GWEN  PROGRAM  AT  NRL  ON  JUNE  7,  1983 

ITS  FUNCTION  IS  TQ  PLOT  THE  HIGHEST  RECORDED  NOISE 
PEAK  FOUND  IN  A  GROUP  OF  10  CONSECUTIVE  DATA  BECORDS 


COMMON /l-EADER/  IDRECOO).  IST(4i.  STSC.  SPSC 
COMMON/DATA/  IDATAt 2, 10000) . NCH 
COMMON/AUTO/  TPD<2. 3. 5000) 

COMMON/TIME/  IH, IM, SEC. ATIMt 5000) 

CHAR AC  TER*30  FILENAM 
CHARACTER  BUFFER *3276S 
BYTE  JDATAI 32748 ). TEST ( 100 ) 

EQUIVALENCE  ( JDATA,  BUFFER ) 

DIMENSION  FMAX<2,  2000) 

DIMENSION  TPDI < 3 ) .  X ( 5000).  PERCENT <2,  5000) . S'GMA(2.  2000) 
DIMENSION  AVE(2.  2000 ).  FT IM<  5000 ) ,  FDATAI2,  5000 ) . TIM (2000) 

LOGICAL  ERRFLG 

OPEN  THE  TAPE  FILE  ON  THE  SCRATCH  DISC  CONTAINING 
THE  DATA  FROM  THE  GWEN  OBSERVING  RUN 

TYPE*.  'INPUT  TAPE  FILE  NAME' 

READ <5. 333)  FILENAM 
FILENAM-'SIAl.  C ODENWALD ]B  DAT' 

FORMAT <5X.  A30) 

0PEN<UNIT-1 ,  CARR  I ACECONTROL* 'LIST ' . 

1 STATUS* ' UNKNOWN ' . RECORDS I ZE*20040  > 

IDATA  CONTAINS  THE  INTEGER  VALUES  FOR  UP  TO 
10.  000  DATA  SAMPLES  FOR  TWO  CHANNELS  EACH 
HAVING  N1  AND  N2  SAMPLES  RESPECTIVELY 

ITPD  IS  THE  VECTOR  CONTAINING  THE  THREE 
THRESHOLD  SENSITIVITY  LEVELS  FOR  THE  TPD 


PROGRAM  BEGINS 

DEFINE  NUMBER  OF  WORDS  IN  DATA  RECORD 

ID  RECORD 

NDA  =  10000 
NID  -  20 
NTOT-NDA+NID 
NBY  »  2*NTQT 

NUMBER  OF  TAPE  RECORDS 


type*,  'type  number  of  records  to  read' 

read*,  mrec 

MREC“1 100 

NSKIP  -  0 

DO  9S  1*1,  30 

IDREC ( I)»0 

CONTINUE 

calibration  constant  for  gain  1  to  convert 

from  data  units  (volts)  squared  to  watts/meter  square 
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M  ui 


pcel-8  t>2#-19 


TYPE*.  ' TAPE  DUMBER ' 

READ*.  I DREC  <  30 > 

;  TYPE*.  'NUMBER  OP  RECORDS  TO  SKIP' 

READ*. NSKIP 

type*,  'NUMBER  OF  RECORDS  TO  (OCESS ' 
read*,  mrec 

TYPE*,  'SKIP  *  '.NSKIP,  PROCESS  =  '.MREC 
IDREC ( 30>«1 

INITIALIZE  PLOT  UTILITY 
CALL  PLOTST  <  5.  'IN'.O) 

;  INITIALIZE  DATA  ARRAYS 

DO  35  1-1.  2 
DO  40  J-I.  100 
PERCENTII,  J)«0. 

DO  41  K*l.  3 
TPD!  I.  K,  J>»0 
*1  CONTINUE 

AC  CONTINUE 

C;  CONTINUE 

DO  76  1*1,  100 
TEST  < I )  *  0 
CONTINUE 

READ  TAPE  RECORD  OF  LENGTH  NBY  BYTES 

ISTT-1 
IPAR  -  0 
NCH  -  2 
I  PLOT- 1 
NNREC-NSKIP 
NR  EC  »  0 

IF (NSKIP  EQ.  0 )  GO  TO  100 
DO  201  II-l, NSKIP 
READ< 1 i S)  BUFFER 
FORMAT (A) 

01  CONTINUE 

MXREC— MREC+NSK IP 
WRITE! 6. 555) 

555  FORMAT!'  TAPE  RECORDS  SKIPPED') 

111  FORMAT! 1H0,  6X,  'TIME'.  6X,  '  RECORD '.  8X,  ' CHANNEL ',  6X, 

1 'AVERAGE '.6X.  ' V-SUB-D/20DB ' ,  10X,  'PEAK  POWER'///) 
IFIRST  *  1 
jgain  -  0 

100  NREC-NREC+1 

IF!NREC. GT.  MREC)  GO  TO  1005 
READ! 1, 5.  END- 1004)  BUFFER 
C  WRITE ! 6, 888)  ! JDATA! I ) .  1-1 ,  10) 

CS88  FORMAT! 1015) 

C  NOW  TEST  THE  FIRST  100  BYTES  TO  MAKE  CERTAIN 

•;  THAT  THE  CURRENT  RECORD  AND  THE  PREVIOUS 

C  RECORDS  ARE  NOT  IDENTICAL  IF  THE  MULTIPLEXER 

I  HANGS  UP  THIS  WILL  BE  THE  CASE 

IF! IFIRST  EQ  0)  GO  TO  889 


53 


THIS  IS  THE  FIRST  RECORD  FROM  THE  TAPE 
DO  87  IM=1.  100 
TEST ( I ) =JDAT A ( I ) 

CONTINUE 
IFIRST  =  0 
GO  TO  887 
DO  88  1=1,  100 

IF ( TEST ( I )  NE.  JDATA(I))  GO  TO  887 
CONTINUE 

IF  ALL  ELEMENTS  NOT  IDENTICAL  THIS  IS  GOOD  DATA 
WRITE!  (a,  886)  NREC 

FORMAT (  '  DETECTED  REPEATED  RECORD  '.14,  '  RECORD  SKIPPED') 

NREC=NREC-1 

GO  TO  100 

CALL  DECODE! JDATA,  NBV,  NREC,  ISTT) 

set  if  the  gain  setting  has  changed 

IF ( JGA IN.  EG  IDREC ( 16 ) )  GO  TO  9G7 
WR ITE (  6,  938 )  IDRECU6) 

FORMAT! '  THE  GAIN  IS  NOW  SET  TO  ',15) 

UCAIN= IDREC! 16) 

COMPUTE  AVERAGE  VALUE  AND  SIGMA  OF 
DATA  IN  RECORD 

PLOT  THE  FIRST  RECORD  IN  CHANNEL  1  AND  2 
IF  THIS  IS  THE  FIRST  RECORD  PROCESSED 
EXPRESS  X  AXIS  AS  TIME  IN  MILLI  SECONDS 

IF! IPLOT  NE.  1 )  GO  TO  50 
DELT-1000  /6200. 

SUM 1=0 
SUM2-0 
KK*1 

DO  32  1=1,  5000 
FTIM! I >=I*DELT 

fdata!  1, i  )  =  ( 1  0*1  data ( 1 ,  i  )  / 16.  >**2. 
f data <2, i )»( 1.  0*IDATA<2,  I > / 16.  )**2. 

SUM1-SUM1+FDATA! 1,  I) 

SUM2»SUM2+FDATA!2,  I) 

CONTINUE 
AVEKK= SUM 1/5000. 

AVELL=SUM2/5000 
XST=10.  /! 5000.  *DELT ) 

YS1  =  2/AVEKK 
YSO  =  2/AVELL 

WRITE! 6. 89)  ! IDREC !KKK),  KKK=1,  20) 

FORMAT!'  TAPE  ID  '  ,10110/) 

IF(NREC.  EG.  1 )  WRITE!6.  Ill) 


determine  the  g  axis  limit  IN  UNITS  OF 
10!-10>  WATTS/METER  SQUARE/HERTZ 

one  channel  is  at  a  gain  of  1  while  the 

second  channel  is  30db  lower  in  gain  (factor  of  33) 

SCALE-YS1 

I F ( F DA T A < 1 ,  10) .  LT.  FDATA! 2,  10) )  SCALE=YSO 
BBB-8  *pc  a  1*1 .  Oe 1 5*< 10** (30.  /20  ) ) /SCALE 
AAA=8  *pcal*l  0*15/SCALE 
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n  o  o 


cl  *  aaa/3100 
c2  -  bbb/3100 

if <fdata<2,  1 >  9 1.  f data (1  -  1 ) >  cl  »  bbb 
t  f !  f da  ta(2.  1)  gt  fdata!l.  1))  c2  «  aaa 
CALL  PLOT <  5.  2.  ,  -3) 

CALL  PLOTTlFTIM,  FDATA,  3.  XST,  SCALE.  5000,  1. cl,  c2> 

IPLOT-O 

DO  150  J-1.2 

CALL  S  T  AT ( J ,  SIC,  AV, AAV,  10000, NREC. FAX  > 

FMAX  <  J, NREC) -FAX 

FDATA(2, NREC)-! 1.  0*IDATA<2,  1000) )  **2 
IvO  CONTINUE 

30  TO  100 
POO  WRITE (6, 2) 

2  FORMAT (  '  ERROR  IN  INITIALIZING  TAPE') 

GO  TO  1000 
*  0 1  WRITEld, 3) 

j  FORMAT (  '  END  OF  TAPE  ASSUMED  '  ) 

GO  TO  1005 

1 004  WRITE! b. 1006)  NREC 

lCOa  FORMAT!'  END  OF  DISK  OR  TAPE  FILE  DETECTED 
1  AFTER  RECORD  NUMBER  '.15) 

1005  NREC  -  NREC-1 

0  NOW  COMPUTE  AVERAGES  FOR  EACH  GROUP  OF  TEN 

0  RECORDS 

YMAX1=0. 

YMAX2-0. 

DO  412  K-l,  2 
J»1 

AT-ATIMU  ) 

SMX-  0 

DO  411  1-1,  NREC 

IF ( ABS! ATIM! I )-AT).  GT.  .  01  )  CO  TO  413 
IF! FMAX <K.  I).  GT.  SMAX5SMX  -  FMAX (K,  I) 

IF ! FMAX ! 1 ,  I )  GT.  YMAX1)  YMAX1-FMAX (1,1) 

IF (FMAX (2.  I ) .  GT.  YMAX2)  YMAX2-FMAX ( 2,  I) 

IF! I  EG. NREC)  GO  TO  413 
GO  TO  411 

-113  FMAX (K. J) -SMX 

FDATA! 1, J) -ATIM! 1-1 ) 

FDATA! 2, U  > -FDATA (2.  I  —  1 ) 

AT-ATIM! I ) 

SMX  ■  FMAX (K, I) 

J-U+l 

•HI  CONTINUE 

J-J-l 

AVE1K.  UJ-SUMl/lO. 

SIGMA1K, J>-20.  *ALOG 1 0 < AVE ( K.  J) /SIGMA1 ) 

FMAX (K. U) -SMX 
FDATA! 1, U)- AT 
412  CONTINUE 

WRITE! 7, 434)  IDREC(5),  IDREC!4),  IDREC(3),  1ST! 1 ),  I ST <2 
434  FORMAT!'  '.514) 

DO  415  1-1 1 J 

:  TYPE*.  I, FMAX! 1,  I ) ,  FMAX !2.  I ),  SIGMA! 1,1),  SIGMA<2.  I  ) 

WRITE! 7, 435)  I. FDATA! 1,  I ),  FMAX (l,  I >.  FMAX (2,  I  ) 

425  FORMAT!'  '.14.  3FI5.  2) 

ATIM! I ) -FDATA! 1 , I ) 

415  CONTINUE 
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XS3  -  10  / 30 
XS3-10  /80 


COMPUTE  Y  AXIS  SCALES  AND  LIMITS 
YS1  -  1  /SMT 1 

IF ( AWE (I.  10 )  LT  AVE (2,  10  > )  YS1  -  1  /SMT2 

BBB-S  *3.  08E-9/YS1 
AAA-8  *2.  02E-4/YS1 
TYPE*,  YS1,  BBS . AAA 
Cl  «  SBB 
C2  *  AAA 

IF(AVE<1,  10)  LT  AUE<2,  101)  C2  «  BBB 
IF < AVE ( 1,  10)  LT  AVE ( 2,  10)1  Cl  «  AAA 
CALL  PLOT ( 11  ,0  .  -3) 

SCALE-4  /2  5*8 
C 1-8  62E-9*4,  /SCALE 
C2-8  62E-9*10**<30.  /20.  )*4.  /SCALE 
IF<  YMAX2  GT  YMAX1)  Cl-8.  &2E-9*4.  /SCALE 
IF(YMAX2.  GT.  YMAX1)  C2-8.  62E-9* 10** ( 1  5>*4.  /SCALE 
CALL  PLOTT ( ATIM,  FMAX,  6,  XS3,  SCALE,  J,  J,  cl,  c2> 

NOW  PLOT  TAPE  ID  PARAMETERS  ON  PLOT  HEADER 

IDREC ( 29) -NREC 
IDREC ( 28)  -  IPAR 
IDREC (27)  -  1 
CALL  PLOT <9  5, 0,-3) 

CALL  HEADER (NID) 

CALL  PLOTND 
CLOSE ( UNI T-l) 

STOP 
pun 


SUBROUTINE  PlOTT' X. V, INDX, XSCALE. YSCALE. NDAT. NREC. chi, ch2) 

THIS  is  the  subroutine  package  which  generates  plots  OF  the 

INDX  »  1  AMPLITUDE  PROBABILITY  distribution 

INDX  «  2  TIME  PROBABILITY  DISTRIBUTION 

INDX  =  3  PLOTS  AMPLITUDE  FOR  NDAT  DATA  SAMPLES 

INDX  =  4  PLOTS  SIGMA  US  RECORD  NUMBER 

INDX  ■  5  PLOTS  AVERAGE  AMPLITUDE  VS  RECORD  NUMBER 

INDX  «  6  PLOTS  MAXIMUM  DEVIATION  VS  RECORD  NUMBER 

WRITTEN  BY  DR  5TEN  ODENWALD 

NAVAL  RESEARCH  LABORATORY 
CODE  4138-0 

WASHINGTON,  D  C  20375 


CHARACTER*10  XTIT. YTIT, XTIT1. YTIT12, YTIT3 
CHARACTER*10  YTIT4,  YTIT5,  XT IT2, YTITfe, XTIT99 
CHAR AC  TER *1  VAR 
CHAR AC  TER *30  XTJT34 
CHARACTER*12  TITLE1,  TITLE2 
COMMON/HEADER/ 1 DREC <301,  IST<4), ST, SP 
COMMON/TIME/ IH,  IM,  SEC,  AT IM< 50001 
c ommon /au t o/ tp ( 2,  3,  5000) 

DIMENSION  X  <  5000 ) ,  Y  ( 2,  5000) 

REAL*4  FI  I 
INR-1 

3333  CALL  PLOT <8.  >0.  , 2 ) 

CALL  PLOT <  9.  ,  10,  ,2) 

CALL  PL0T(0.  . 10  ,  2) 

CALL  PLOT  ( 0.  ,0.2) 

IF<INDX.NE.  UCALL  PL0T14  ,  0.  .2) 

IF( INDX. NE  1 )  CALL  PL0T<4.  .  10.  . 2) 

CALL  PLOT  (  0.  ,0,3) 

0  INITIALIZE  TITLES 

XTIT1-  '  POWER 
XTIT2-  •  TIME 
XTIT34* '  ' 

YTIT12* '  PERCENT  ' 

YTIT3  »'  POWER 
YTIT4  »'  VD 
YTIT5  *' 

TITLE1* 'FIRST  RECORD' 

TITLE2* 'ENTIRE  TAPE' 

YTIT6»  'PEAK  POWER  1 

5  IF< INDX. NE.  1 )  CO  TO  10 

15  XTIT-XTIT1 

YTIT-YTIT12 
GO  TO  50 

10  IF ( INDX.  NE.  2)  CO  TO  20 

XTIT-XTIT2 
YTIT-YTIT12 
CO  TO  50 

20  IF ( INDX.  NE.  3  GO  TO  30 

XTIT-XTIT34 
YTIT-YTIT3 
GO  TO  50 

30  IF (INDX  NE  4)  GO  TO  40 


o  n  c» 


YTIT-/TIT4 
50  TO  50 

■l'‘,  IFUNDX.NE  5 GO  TO  45 

XTIT-XTIT34 
YTIT-YTIT3 
GO  TO  50 

45  IF< INDX  NE.  6)  GO  TO  46 

XTIT-XTIT34 
YTIT-YTIT6 
GO  TO  50 

4e  WRITE (6.1) 

1  FORMAT  <  '  NO  CODE  FOR  THIS  PLOT.  SELECT  ONE') 

READ( 5.2)  INDX 

C  FORMAT (II) 

IF< INDX  EQ  0)  RETURN 
GO  TO  5 

:  NOW  DRAW  TIC  MARKS  QN  X  AXIS  AND  LABEL 

50  IF< INDX  EQ.  3)  GO  TO  60 

1F< INDX. EQ.  4)  GO  TO  60 

IF ( INDX  EQ.  5)  GO  TO  60 

IF ( INDX  EQ. 6)  GO  TO  60 

C  DRAW  LOGARITHMIC  INTERVALS 

:  ON  X  AXIS 

: 3  DO  100  1-1. 3 

DO  101  J-l.9 
FII-J*10**(I-1) 

X 1-XSC ALE*ALOC 10  < F 1 1 ) 

CALL  PLOT <8.  . XI.  3) 

CALL  PLOT <  7  9,  XI,  2) 

IF(J  EQ.  t>  CALL  PLOT ( 7  S.  XI. 2) 

.Cl  CONTINUE 

.00  CONTINUE 

C  NOW  LABEL  DECADES 

\- 

DO  150  1-2.4 
X2-XSCALE*< 1-2) 

IF  < INR  EQ.  1 ) CAL  I B— CHI 
IFUNR  EQ.  2)  CALIB-CH2 
IF< INDX. EQ  2 )  CALIB-CH1 
XN-CALIB*10**( 1-2) 

CALL  NUMBERO  2.  X2.  l.XN.  90  .5) 

150  CONTINUE 

CALL  SYMBOL (9  ,  3  5.  2.  %REF ( XTIT > .  90.  ,  10 ) 

IF < INDX. NE.  3)  GO  TO  70 

LABEL  TVC  PLOT  BY  THE  RECORD  NUMBER  RETRIEVED 
XTIT- 'RECORD' 

CALL  SYMBOL <9  ,  3.  ,  .  3.  7.REF  (  XTIT ) .  90.  .  6 ) 

CALL  NUMBER  <  9  .5  5.  .  3,  FLOAT ( NREC ) .  90.  .  1 ) 

GO  TO  70 

DRAW  X  AXIS  FOR  LINEAR  TIME  PLOT 

IF( INDX  NE  3)  GO  TO  61 
00  200  1-1.9 
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Xt-I*l 

DEC  *  X  < 1 ) *500 
X2  -  D£C*X 1 
CALL  PLOT (8  .  XI,  3) 

CALL  PLOT <7  9,  XI, 2) 

CALL  PLOT  <  4  ,  XI,  3) 

call  plot < 3,  9, xi. 2) 

CALL  NUMBER (8  2,  XI-  2,  1,X2,  90  -1) 

00  CONTINUE 

30  TO  63 

i  DEC-30  /tO 

mark  th#  x  axis  in  hours 

DO  202  1-1, 9 
X  1  - 1  *  1 

12  -  1ST i 1 >-rDEC*Xl 
IF ( 12  ST  24)  12-12-24 
CALL  PLOT ( B  ,  XI,  3) 

CALL  P LOT <  7  9,  XI,  2) 

CALL  PLOT (4.  ,  XI, 3) 

CALL  PLOTO.  9,  XI,  2) 

CALL  NUMBER (8.  2,  XI-  2,  .  1, FLOAT! 12) ,  90.  .  1 ) 

02  CONTINUE 

3  XTIT- ' TIME  IN- 

CALL  SYMBOL <9.  0,  3  ,  2,  ‘/.REF!  XTIT).  90.  ,  7) 

IF< INOX. NE  3)  00  TO  201 
XTIT- 'MILL I  SECO' 

CALL  SYMBOL (9  0,  4  6,  2,  */.REF!  XTIT) ,  90  .  10) 

XTIT-'NDS' 

CALL  SYMBOL  (9  0,  6.  6,  2,  */.R£F<  XTIT) ,  90  ,  3) 

01  XTIT- 'RECORD' 

IF( INDX  NE.  3)  CALL  SYMBOL!-.  2.  3.  .  .  3,  %REF(TITLE2> . 90.  ,12) 
IF!  INDX.  EG.  3) CALL  SYMBOL!-  2.  3.  .  .  3.  7.REF!XTIT),  90  .  6) 
XTIT- 'HOURS  ' 

1FUNDX.NE.  31CALL  SYMBOL! 9.  0.  4  6.  .  2,  V.REFUTIT).  90.  .  7) 

IF< INDX. EQ. 3)CALL  NUMBER!-  2.  5,  5,  .  3.  FLOAT  ( NREC  > .  90.  .  1  > 


NOW  LABEL  Y  AXIS 

0  IF< INDX.  EQ.  3)  00  TO  80 

IF! INDX.  EQ.  4)  00  TO  80 
IF!  INDX.  EQ.  S)  00  TO  80 
IF! INOX  EQ.  6)  CO  TO  80 

DRAW  LOO  TIC  MARKS  ON  Y  AXIS 

DO  250  1-2.6 
Y2-8.  -YSCALE* ( 1-1 ) 

YDX-t-3 

IF! INDX.  EQ.  1)  CO  TO  239 
IF!  I. LT. 3)YDX-I-2 
59  YN«10**YDX 

CALL  NUMBER !Y2. -  5.  1. YN, 90. ,3) 

50  CONTINUE 

DO  300  1*1.  6 
YO-8.  -YSCALE*! 1-1 ) 

DO  301  J-l. 10 
R2-J 

/ 1 -YO- YSC ALE*ALOC 1 0 ! R2 ) 
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CAUL  r  LOT (VI  .0  .3) 

CALL  P  LOT  <  Y 1 .  1.2) 

CONTINUE 
CONTINUE 
GO  TO  90 

DRAW  LINEAR  r  AXIS  SCALE 

FMAX  -  4  /yscal* 
if(ind«  »n  3)  findi  -  chi 
ifUndi  »q  5;  £max  *  chi 
IF < INDX.  EQ.  6)  FMAX"  CHI 
DO  400  1-1. 20 
Y 1  —  4*1 

CALL  PLOT ( Yl, 0  , 3) 

CALL  PLOT< Yl,  1. 2) 

CONTINUE 

LABEL  THE  Y  AXIS 

CALL  NUMBER  ( 0  ,  -.  7,  .  1.  FMAX,  90.  ,  5) 
CALL  NUMBER <3.  9,  -  7.  .  1,  0.  .  90.  ,  5) 

IF < INOX. EQ.  3)  FMAX  -  ch2 
IFUNDX.  EQ.  5>  FMAX  -  ch2 
IF < INDX.  EQ.  6)  FMAX-CH2 
CALL  NUMBER (4  1.  -  7,  .  1,  FMAX,  90.  ,  5) 
CALL  NUMBER < 9.  ,  -  7,  1.  0.  .  90.  .  5) 

LABEL  CHANNELS  ON  Y  AXIS 


CALL  SYMBOL  ( 4  5.  -1.  2,  .  2,  XREF<  YTIT) ,  -180.  ,  10) 

ROT— 180. 

51- 4.  5 

52- 3.  6 

53- 4.  3 

54- 3.  0 
S3—.  9 
S6— 1. 

IF< INDX  EQ. 2)  GO  TO  920 
IF ( INDX. EQ.  5)  GO  TO  910 
IF ( INDX.  EQ  3)  GO  TO  910 
IFIINDX  EQ.  1)  R0T*90. 

YTIT3-U0  WATTS' 

CALL  SYMBOL (SI ,  S5.  .  I,  XREF (YTIT3 ) ,  ROT,  10) 

IFUNDX.  EQ.  1)  CALL  SYMBOL <9.  2,  3.  4,  .  I,  XREF <  YTIT3 ) ,  ROT ,10) 
YTIT4-  '/METER ' 

CALL  SYMBOL (S2,  S5,  .  1, XREF (YTIT4) , ROT,  6) 

IFUNDX.  EQ.  I  )CALL  SYMBOL  (9.  2,  4.  4,  .  1 , 7.REF  ( YTIT4  ) ,  ROT,  6) 
YTIT* ' -15 ' 

IF< INDX.  EQ.  6)  YTIT* ' -10 ' 

CALL  SYMBOL (S3.  S6.  .  1.  XREF(YTIT), ROT, 3) 

IF<  INDX.  EQ.  I  JCALL  SYMBOL (9.  1,  3.  6,  .  1,  '/.REF<  YTIT) .  ROT,  3) 
YTIT" ' 2 ' 

CALL  SYMBOL < S4,  S6,  .  1 , XREF (YTIT) ,  ROT.  1) 

IFUNDX.  EQ.  DCALL  SYMBOL  <  9.  1 ,  5.  .  1 , '/.REF  (  YT  IT  ) ,  ROT,  1  ) 
IFUNDX  EQ.  1)  GO  TO  93 
YTIT- 'CHI ' 

CALL  SYMBOL (2  ,  -l.  2,  .  2, XREF(YTIT), -180  .  3) 

YTIT- 'CH2 ' 

CALL  SYMBOL ! 6  , -l. 2,  2, XREF(YTIT), -1Q0.  , 3) 


NOW  PLOT  POINTS 


IF <  INDX  EQ  3)  GO  TO  600 
IF ( INOX  EQ.  4)  GO  TO  600 
IF < INDX  EG.  5)  GO  TO  600 
IF < INDX  EG  6  >  GO  TO  600 

C  SET  UP  FOR  LOG  PLOT 

C  FIRST  CHECK  TO  SEE  IF  THE  TPD  IS  TO  QE  PLOTTED 

IF < IND X , NE.  2)  GO  TO  95 
DO  97  12-1,  3 
YO  *  4 

DO  96  1 1-1. 2 

IF <  1 1  EG.  2)  YO  =  8 

A— TP (II, 12, 1) 

IF  <  A.  LT  .  001)  A  -  001 

XX-O.  0 

IF  <  X  < l > . GT.  0.  >XX  -  XSCALE<<AL0G10<  X ( 1 ) >  — 1.  ) 

YY  -  Y0-YSCALE<<3<AL0G10(A> > 

CALL  P LOT <  YY ,  XX,  3> 

DO  98  1*2,  NDAT 
A-TPUl, 12, I) 

IF < A.  LT.  .  001)A  -  .  001 
XX-XSC ALE<<ALQG10<  X  < I ) ) -1 .  ) 

YY-YO-YSCALE< < 3<Al0G 10 < A ) ) 

CALL  PLOT <  YY,  XX, 2) 

IF< I.  EG. 2 ) CALL  SYMBOL< YY. XX,  .  2,  12,  90.  , -1 ) 

98  CONTINUE 

96  CONTINUE 

97  CONTINUE 
CO  TO  94 

95  C0£-AL0C10<  X  < 1 ) ) 

XX-XSC ALE<<AL0010<  X  < 1 ) >-COE> 

A-YUNR.  1) 

IF<A.  LT  000001)  A  -  000001 

YY-8.  Q-YSCALE* <6<AL0C10<  A) ) 

CALL  P LOT < YY , XX,  3) 

DO  650  1-2. NDAT 
A  -  Y<  I Ml.  I) 

IF  <  A.  LT.  .  000001 >  A  -  000001 

XX-XSC ALE<<AL0G10<  X  < I ) J-COE) 

YY-8.  0- YSCALE<<6<AL0G10< A) ) 

CALL  P LOT < YY , XX, 2) 

630  CONTINUE 

94  CALL  PLOT <0.  .  0.  .  3) 

IF < INDX. EQ.  2)  RETURN 
XTIT99- 'CHANNEL  ' 

CALL  SYMBOL <-.  2.  4.  .  2,  XREF <  XTIT99) .  90.  ,9) 

CALL  NUMBER < - .  2.  6.  .  2,  FLOATUNR) ,  90.  ,  1> 

I F < I NR .  EG.  llCALL  PLOT< 11.  . 0.  .  -3) 

IF < INR. EQ.  2)  RETURN 

INR-INR<1 

GO  TO  3333 

600  XX-XSC ALE*X  < 1 ) 

YY-4. -YSCALE<Y< 1,  1 ) 

IF<YY.  LT.  0.  )  YY  ■  0. 
call  p lot  <  yy,  ii,3) 

IF < INDX.  EG.  6 >CALL  PLOT < 4.  , XX, 3 > 

CALL  PLOT< YY,  XX,  2) 
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•  •. .  ,  i  <j.  ,\av,  r-OAi .  . o 

IMIS  Rnur  I  lE  TAKES  THE  DATA  FROM  EACH 
OF  TWO  DATA  CHANNELS  AMD  COMMUTES  THE 

AVERAGE  value  and  sigma  for  a  single  data 
RECORO  CONSISTING  OF  5000  POINTS  PER  CHANNEL 

WRITTEN  3Y  DR.  STEN  ODENWALD 


COMMON/TIME/  IH.  IM,  SEC.  TIM< 5000) 

COMMON /DATA/ I  DATA (2.  10000) . NCH 
common/ header/ idrec (30).  itd(4),  sts,  sps 
DIMENSION  FDAT(2, 10000) 

I  COUNT  -0 

SUM-0 

SUM2-0 

C 

c  read  the  gain  setting  and  compensate 

c  idrec(16>  «  1  means  0  db  so  gain  -  t. 

c  idrec  <  16)  «  2  mean*  24  db  so  gain  -  .0631 

c  idrec  (16)  »  4  means  43.5  db  so  gain  -  .00668 

c  idrec! 16)  *  8  means  63  db  so  gain  -  .  000708 

c  if  idrec(16>  equals  anything  else,  gain  *  1. 

c 

igain-IOPEC! 16) 
gain  -  1. 

ifdgain.  eq.  2)  gain  *  .0631 
i  f  <  iga  i  n.  eq.  4)  gain  -  .00668 
i  f  (  iga  in.  eq.  8)  gain  «  .000708 
C  DO  90  1*1.  5000 

I  FOAT  <U,  I >»<  < 1.  Q*IDATA{ J.  I)/<16.  )>**2.  )/GAIN**2 

'CO  CONTINUE 

COMPUTE  MEAN  SQUARE  VOLTAGE  -  AVERAGE  POWER 

:  DO  100  1-1.5000 

:  SUM-SUM+FDAT! J.  I ) 

:  SUM2="SUM2+SQRT<FDAT<  J.  I)) 

IlOO  CONTINUE 

C  AV— SUM / 5000. 

C  AAV-SUM2/5000. 

:  TYPE*.  AV 

C  COMPUTE  SIGMA 

C  SIGMA-O. 

C  DO  200  JJ-1.  5000 

:  SIGMA-SIGMA+<AV-FDAT< J.  JJ)  >*< AV-FDAT! J,  JJ) > 

:200  CONTINUE 

C  SIG-O. 

C  IF(ABS<SIGMA).  LT.  .  0001)  RETURN 

C  SIG-20. «L0G10< SORT ( AV) /AAV) 

FMAX  »  0. 

DO  250  JJ- 1.5000 
C  FM-FDATf J,  JJ) /av 

if  (idrec<4).  ne.  6)  go  to  251 
*  if(idrec<5>.  It.  15)  go  to  251 

f .n-  Id  i  ta<  j.  j  j  ) 

CM  IF  (FMAX.  LT.  FM)  GO  TO  201 

CO  TO  250 
^  .  1  Fhiax  «  FM 

l-C  CONTINUE 

r  MAX »FMAX*FMAX / (256.  *GAI N ) 

WR  I TE (  6.  10)  IH.  IM.  SEC.  NREC.  J.  AV,  SIG,  F"AX 
10  FORMAT  (  IH  .  12.  *  ',12,'  ' ,  F6  3.  1 10,  5X.  1 10. 

15X.E10  2.5X.F10  2.  5X,  1E10  3) 

RETURN 

END 
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SUBROUTINE  AUTQi  ITPD.  X,  TPDi NX, NRC > 


AUTOCORRELATES  the  ARRAY  ITPD  OF  LENGTH  NRC 
WRITTEN  BY  DR  STEN  ODENWALD 


COMMON /DATA/ I DATA (2,  10000 > . NCH 
DIMENSION  TPD<2, 3. 5000), X<5000>, ITPD13) 

NORMALIZE  THE  SIGNAL  STRENGTH  TO  UNITY 
AMPLITUDE  FOR  EACH  POINT  ABOVE  THE 
THRESHOLD  ITPD  AND  COMPUTE  AUTOCOR 

A  =  5000* < NRC- 1 ) 

B  -  5000*NRC 
DO  700  L* 1 , 3 
DO  400  J- 1 , 2 
DO  600  1-1, NX 
TD  -  0. 

IMAX  «  5000- X ( I ) + 1 
DO  550  K«l,  IMAX 
INDX*X  <  D+K-l 

IF< IDATAI J,  K)  .  LT  ITPD(L))  GO  TO  550 
IF< IDATAIJ,  INDX)  LT  ITPDIL) >  GO  TO  550 
TD  »  TD+1. 

CONTINUE 

TPDI J. L. I )  «  (TPDIJ, L,  I)*A+TD)/B 

CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

END 
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SUBROUTINE  DECODE ( I ARRAY , NBVTE, NREC, 1ST) 


THIS  PROGRAM  UNPACKS  THE  DATA  UALUE  FROM 
THE  2-BYTE  WORD  REPRESENTATION  AND  STORES 
THE  RESULT  IN  'ARRAY' 

WRITTEN  BY  DR-  STEN  QDENWALD 


COMMON/TIME/  IH, IM, SEC, TIM< 5000) 
COMMON/HEADER/ IDREC< 30 )  ,  ITD<4>. STS, SPS 
COMMON/DATA/IDATAI 2,  10000) 

INTEGER*2  I ARRAY  < 10050  > 

INTEG-N0YTE/2 

SWAP  THE  BYTES 

CALL  BYSWAP ( IARRAY.  INTEG) 

1 1-0 

DO  100  1=1.10000,2 
II-II+l 

I DATA ( 1,11 l-IARRAY ( I  ) 

CONTINUE 

II-l 

DO  ISO  1-2.10000.2 
I DATA (2, I I )«I ARRAY ( I ) 

II-II+l 
CONTINUE 
K— 1 

DO  200  1-10001. INTEG 
IDREC  <  K )«I ARRAY  < I > 

K-K+l 

CONTINUE 

CALCULATE  TIMES 


A- IDREC ( 1 )  *.  t 

IF< IDREC< 1 ) .  LT.  0)  A-3276.  8+< 32768+IDREC < 1 > > *  1 
B-2.  + 1 DREC ( 2 ) *3276 .  8 

TS»( A) +IDREC ( 8  J  +  IDREC  <  7  > *60.  + 1 DREC (6) *3600.  +B 

IH-TS/3600 

FH-TS/3600. 

IM«60*(FH-IH) 

FM-60.  *<FH-IH> 

SEC-60. *<FM— IM) 

IF( 1ST.  PC.  1)  GO  TO  120 

ITD< 1 )-IH 

ITD<2)-IM 

STS-SEC 

I ST-0 

ITD<  3 ) -IH 

ITD<4)-IM 

SPS- SEC 

TIME  IN  HOURS  SINCE  THE  START  OF  THE  TAPE  IN  MIN  PAST 
START  HOUR 

TlM<NREC)-( IH-ITD( 1 ) )+< IM+(SEC/60.  ) >/60. 

RETURN 


SUBROUTINE  HEADERtN) 


PRINTS  A  HEADER  AT  END  OF  PLOT 
WRITTEN  BV  DR.  STEN  ODENWALD 


COMMON / HEADER / I DREC  (  30 )  ,  I ST ( 4 , ,  STSC .  SPSC 
CHAR AC  TER *40  XTIT ( 30 ) .  IMO ( 12 > . COL 
XTIT < 1 >  = 'SUMMARY  OF  GWEN  DATA  TAPE' 

*  tit (2)  =  ' FOR ' 

* t i t ( 3 ) = 'TAPE  NUMBER ' 

* t i t ( 4 ) ■'MONITORING  FREQUENCY ' 

x  t i t  <  5 ) * ' C  HANNEL  l' 

x  t  i  t  ( 6  >  * '  C HANNEL  2' 

x  t  i  t  ( 7  > « 'BANDWIDTH ' 

x t i t  <8  >*'HZ ' 

x  t i t ( 9 )  =  'LOCATION ' 

x t i t < 1 0 ) »  ' MAR YL AND  POINT' 

xtittll>='START  TIME' 

xtitt 12)« 'STOP  TIME' 

xtit( 13)- 'NUMBER  OF  RECORDS' 

xtitU4>-' 

x  t  i  t  <  1 5  )= ' 

x t i t < 16  >  = 'SAMPLE  DURATION' 

xtitt 17)-' INTERVAL  BETWEEN  SAMPLES' 

XTIT(18)«'SEC' 

imotl)-'  JANUARY  ' 

imo(2>«'  FEBRUARY  ' 

imo  <  3) *  '  MARCH 

into  (  4  )  =*  '  APRIL  ' 

imo<5)»'  MAY 

imo(6) * '  JUNE 

imo  (  7 )  *  '  JULY 

imo(8>«'  AUGUST 

imo  <  9 ) *  '  SEPTEMBER' 

imotlO)*'  OCTOBER  ' 

imo(U)*'  NOVEMBER  ' 

lino (12)-'  DECEMBER  ' 

col*':  ' 

call  p  lot<  1.  0.  0.  .  -3) 
call  p  1  ot  ( 9.  ,  0.  ,  2) 
call  p lot < 9  .  10.  . 2) 
call  p  lot(0. . 10. . 2) 
call  p  lot  (0.  .0.  .2) 

call  symboltl  .  1 .  .  .  3.  Xr  ef  <  x  t  i  t  <  1 )  )  <  90.  ,  25) 
call  symbol  <2.  .  2.  5.  2.  7.r«F<  x tit (2>  ).  90.  .  3> 
call  numbar  <2.  >4.  .  .  2<  f  lost  <  idree  ( 5 )  >  >  90.  >  1> 
call  symbol<2.  >.S.  2.  .  2<Xr»f<imo(idr*c(4))),90 
call  numbar <2.  <7  Si.  2ifloat(idrac<3>)i90.  •  1) 
call  symbol  (3.  •  1  5.  .  3.  -/.raf  (  xti  t  <3> ) ,  90.  .  1 1 ) 
call  numbar<3.>7  .  .  Si  f 1  oat < i dr ac (30 ) ) ( 90.  ■  1 ) 
call  symbol  <4.  ,  1  5.  .  2.  */.r af  (  x  t i  t  ( 4 )  ) .  90.  .  21 ) 
ID9»IDREC ( 9) / 100 

call  numbar<4  ,  6.  .  3i  f 1  oat < i d9 > ■ 90.  .  1 ) 
call  s  ymb  o  1  (  4  ,  8  S.  2.  Xraf  <  'KHZ  90.  ,  3 ) 
call  symbol  <4  5.  1.  5.  .  2.  Zraf  (  x  t  i  t  <  7  >> .  90.  .9) 
call  numbar  <  4.  5,  6  ,  .  3i  float(idrac(lO))i 90.  •  1 
CALL  SYMB0L<4.  5,  8.  5.  .  2.  '/.REF <  XTIT ( 8 >  ) .  90.  .  2) 
call  symbol(5.  ■  1.  5i  .  2i  Xr  af  C  xt  i  t  (9)  )  •  90.  iS) 
call  symbol <5  .4  ,  3.  XraF ( x ti t ( 10 ) ) . 90.  .  14 ) 


call  s  y  mb  o  l  (  6  ,1  5.  2,  7.r  «f  <  *  1 1 1 1  1  l  )  ) ,  90  ,101 

IF  (  1ST  (!)  <5T  24)  IST<1)  •  IST(l)-24 

CALL  NUMBER (a.  .4  ,  .  2. FLOAT ( ISTl 1 >  >. 90  -  0) 

CALL  SYMBOL (6  ,5  ,  2,  7.R£F<C0L>.  90  .  1  ) 

CALL  NUMBER  1 6  6  .  2, FLOAT < 1ST ( 2 >  >  <  90  ,0) 

CALL  SYMBOL  ( 6  .7  ,  2, '/.REF ( COL  >,  90  .1) 

CALL  NUMBER(6  .8  ,  2,  STSC,  90.  ,  0) 

call  symbo 1  <  6  5,1  5,  2.  ‘4ref  <  x  t  i  t  (  12  >  > ,  90  ,10) 

IF<ISTi3)  GT  24)  IST(3)«IST<3)-24 

CALL  NUMBER (6  5.  4  .  2, FLOAT < 1ST < 3 ) ) .  90  ,  0 ) 

CALL  SYMBOL  <  6  5.  5  .  ,  2.  '/.REF  (COL  > ,  90  .  ,  1  ) 

CALL  NUMBER  (e>  5,  6.  .  2,  FL0AT(IST(4)  > ,  90.  ,  0) 

CALL  SYMBOL  (6  5-7  .  .  2.  ‘/.REF < COL) ,  90  ,  1  ) 

CALL  NUMBER ( 6.  5.  8  .  2- SPSC. 90  ,  0) 
call  symbo  1  <  7  0.  1 .  5.  .  2.  7.ref  <  x  t  i  t  <  13)  > .  90.  ,  17  > 
call  number (7.  0. 7.  0.  2.  f 1  oat ( l dric < 29 ) ) , 90  .,  1 ) 
call  s  ymb  o  1  <7  5-  1  5,  2.  7.r  *4  (  x  t  i  t  (  1 4 )),  90  -13) 
i f < i dr ec < 28 >  ne  0 )  go  to  10 

call  symbol<7  5- 7.  5<  .  3.  Xref  <  x  t  i  t  <  15) ) ,  90,  .  4 ) 
go  to  20 

call  number  (7  5-  7.  5-  .  2-  f  1  oat  <  l  dr  ec  (28 )  ) .  90.  -  1  ) 

call  symbol  (8.  0.  1.  5.  .  2.  -/.r e # <  x  t i t  ( 16 )  1. 90.  ,  15) 

F-5000. /(6200  ) 

call  number  <8.  0.  7.  5.  2.  F.  90.  -  3) 

CALL  SYMBOL  ( 8.  0.  9  ,  .  2, ‘/.REF<  XTIT  <  18  >  > .  90.  ,  3  > 
call  symbol  (8.  5-  15,.  2,  */.ref  <  x  t  i  t  <  17  >  ) ,  90  ,24) 
F«IDR£C  < 13 ) *.  1 

call  number  (8.  5.75,  2, F,  90.  .1) 

CALL  SYMBOL ( 8.  5.9  .  .  2.  XREF(XTIT( 18) )  -  90.  .  3) 

return 

end 
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SUBROUTINE  STATU.  SIG,  AV,  AAV.  NDAT .  NREC .  FMAX  ) 

THIS  SUBROUTINE  TAKES  THE  DATA  FROM  EACH 
OF  TWO  DATA  CHANNELS  AND  COMPUTES  THE 
AVERAGE  VALUE  AND  SIGMA  FOR  A  SINGLE  DATA 
RECORD  CONSISTING  OF  5000  POINTS  PER  CHANNEL 

WRITTEN  BY  DR  STEN  ODENWALD 


COMMON /TIME/  IH,  IM,  SEC.  TIM( 5000) 

COMMON/DATA/ !DATA< 2.  10000), NCH 
c ommon /header / i d rtc ( 30 ) . itd<4), sts. sps 
DIMENSION  FDAT (2. 10000) 

I COUNT -0 

SUM-0 

SUM2-0 

read  the  gain  setting  and  compensate 
idrec(l6)  -  1  means  0  db  so  gain  -  1. 
idrec<16)  -  2  means  24  db  so  gain  -  0631 

idrec(16>  «  4  means  43.5  db  so  gain  -  .00668 

idrec(16>  -  S  means  63  db  so  gain  -  .000708 

if  idrec(16)  equals  anything  else,  gain  -  1 

igain-IDR£C< 16) 
gain  -  1. 

if  (  igain.  eg.  2)  gain  -  .0631 
if < igain.  eg.  4)  gain  *  00668 

if < igain.  eg.  8)  gain  =  000708 

DO  90  1-1.5000 

FDAT <  J.  I)-<(1.  0*IDATA< J,  I)/<16.  > ) e*2.  >/GAIN 
CONTINUE 

COMPUTE  MEAN  SQUARE  VOLTAGE  -  AVERAGE  POWER 

DO  100  1*1.5000 
SUM— SUM+FDAT  U.  I  ) 

SUM2-SUM2+SQRT ( FDAT  (  J,  I) ) 

CONTINUE 
A V— SUM / 5000 . 

AAV-SUM2/5000 
TYPE*. AV 

COMPUTE  SIGMA 

SIGMA-O. 

DO  200  JU- 1,5000 

SIGMA-SIGMA* ( AV-FDAT <  J.  JJ> ) * < AV-FDAT (  J.  JJ>  ) 

CONTINUE 

SIO-O. 

IF< ABS (SIGMA).  LT.  0001 )  RETURN 
SIG-20. *L0G10<SQRT( AV) /AAV) 

FMAX  -  0. 

DO  250  JJ- 1,5000 
FM-FDATU,  JJ)/av 
if(idrec(4>  ne.  6)  go  to  251 
if < idrec(5).  It.  15)  go  to  251 
fm-f da  t ( j.  j  j ) 

IF<FMAX.LT  FM>  GO  TO  201 
GO  TO  250 
FMAX  »  FM 
CONTINUE 

WRITE! 6.  10)  IH,  IM,  SEC.  NREC. V. AV, SIG, FMAX 
FORMAT*  IH  ,12,  '  '.12,  V  F6  3,  1 10,  5X,  1 10. 

S<.  E10  2.  5X.  F10  2,  5X,  1E10  3) 

RETURN 


o  r>  o  ri  r>  i  t  ci  (')  r»  n  o  ri  r»  c>  r .  u  r»  o  o  o 


PROGRAM  apd 


& 


THIS  PROGRAM  WAS  WRITTEN  BY  DR.  STEN  ODENWALD  FOR 
Th£  GWEN  PROGRAM  AT  NRL  ON  MAV  9.  1983. 

-UNCTION:  TO  READ  A  9-TRACK  TAPE  AND  CALCULATE 

THE  AMPLITUDE  PROBABILITY  DISTRIBUTION  FOR  USER-SELECTED 

SAMPLES  OF  THE  DATA 


COMMON/^ADER/  I DREC  (  30  )  .  IST(4)  ,  STSC.  SPSC 
COMMON/DATA/  I  DATA < 2.  10000 > , NCH 
COMMON/AUTO/  TPD(2. 3, 5000) 

IOMMON/TIME/  IH,  IM,  SEC, ATIM< 5000) 

Char AC TER *30  FILENAM 
CHARACTER*!*  TITLE 
CHARACTER  BUFFER *32768 
BYTE  UDATAC 32768).  TEST ( 100) 

EGUI VALENCE  (JDATA. BUFFER) 

DIMENSION  FMAX(2.  2000) 

DIMENSION  TPDI <3 ) .  X<5000>. PERCENT 1 2. 5000). SIGMA <2.  2000) 
DIMENSION  AVE<2,  2000) .  FTIMI 5000) .  FDATAI2.  5000) .  TIM(2000) 
LOGICAL  ERRFLC 

OPEN  THE  TAPE  FILE  ON  THE  SCRATCH  DISC  CONTAINING 
THE  DATA  FROM  THE  GWEN  OBSERVING  RUN 

TYPE*.  'INPUT  TAPE  FILE  NAME' 

READ<  S.  333)  FILENAM 
FILENAM-'SIAl:  CODENWALDIB. DAT' 

33  FORMAT < 5X.  A30) 

0PEN<UNIT«1.  CARR  I AGECONTRQL- 'LIST '. 

1 STATUS- ' UNKNOWN ' .  RECORDS  I ZE-20040 ) 

INITIALIZE  PLOT  UTILITY 
CALL  PLQTSTIS.  'IN'.O) 

IFIRST-1 


1 


IDATA  CONTAINS  THE  INTEGER  VALUES  FOR  UP  TO 
10.000  DATA  SAMPLES  FOR  TWO  CHANNELS  EACH 
HAVING  N1  AND  N2  SAMPLES  RESPECTIVELY. 

ITPD  IS  THE  VECTOR  CONTAINING  THE  THREE 
THRESHOLD  SENSITIVITY  LEVELS  FOR  THE  TPD 


PROGRAM  BEGINS 


DEFINE  NUMBER  OF  WORDS  IN  DATA  RECORD 

ID  RECORD 

NDA  -  10000 
NID  »  20 
NTOT-NDA+NID 
NBY  «  2*NT0T 

NUMBER  OF  TAPE  RECORDS 

MREC  »  1120 
NSKIP  -  0 
DO  38  1-1 .  30 
I DREC < I )-0 


V  V 


Vv\ 
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r>  p.i  io  u» 
O  O 


-8  CONTINUE 

DO  23  1-1.2 
DO  24  J-l. 1000 
°ERCENT< I. J ) -0 
CJ  CONTINUE 

10  CONTINUE 

TYPE#,  'TAPE  NUMBER 

READ*.  IDREC < 30 ) 

.00  TYPE*.  'WHAT  RECORD  DO  YOU  WANT  TO  APD?' 

READ*,  IREC 

IF< IREC  EQ  0)  GO  TO  1005 
IF ( IREC. GT  MREC )  GO  TO  100 
NSKIP-  IREC-1 
MREC-1 
IDREC ( 30) -1 

READ  TAPE  RECORD  OF  LENGTH  NBY  BYTES 

ISTT-l 
IPAR  -  0 
NCH  -  2 
I PLOT- 1 

IF (NSKIP.  EQ. 0)  GO  TO  202 
DO  201  II— 1. NSKIP 
READ <1.5)  BUFFER 
FORMAT < A) 

1  CONTINUE 

2  READ < 1. 5.  END- 1004)  BUFFER 

WR 1TE( 6. 886 )  (JDATA(I),  1-1,  10) 

CS88  FORMAT (1015) 

ee7  CALL  DECODE! JDATA,  NBY.  IREC.  ISTT) 


in  if  th*  gain  setting  has  changed 


c 

c 


C 


GAIN-1. 

JGA IN- IDREC  (16) 

IF <  JCAIN.  EQ.  2)  GAIN-  0631 
IF < JGA IN.  EQ.  4)  GAIN  -  00668 

IF< JGAIN.  EQ.  8)  GAIN  -  000708 

DELT-1000.  /6200 
SUM 1-0. 

SUM2-0. 

KK-1 

DO  32  1-1.5000 
FTIM< I )-I*OELT 

f data!  1. i )-< 1.  0*idata<  1,  i )/16.  >**2. 
f data ( 2. i )-< 1.  0*IDATA<2,  I )/16.  )**2. 

SUM 1 -SUM1 +FDAT A  < 1 .  I) 

SUM2-SUM2+FDATA  <  2.  I) 

CONTINUE 

AVEKK-SUM1/5000. 

AVELL-SUM2/5000. 

XST-10.  /<  5000. *DELT) 

YS1— . 1/AVEKK 
YSO  -  .  1/AVELL 

determine  the  g  ails  limit  IN  UNITS  OF 
10<-15)  WATTS/METER  SQUARE/HERTZ 
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o  r*  o 


SC  ALE-  *S1 

.F  i FDaTA! l .  10  *  lT  FDATA(2,  10)  )  SCALE-vSO 
jBB-S  »0  62«-4«lCi**!30  /20  > / ( GAIN*SCALE > 
aaa»8  *9  62a-4/ ( Gain*SCALE ) 
cl  »  aaa/3100 
c  2  «  b  0  b / 3 1 00 

if(fdata(2.  1  )  g t  f d a t a ( 1 ,  t  )  >  cl  -  bbb 

1 1  (  f  da  ta  <  2.  1  )  gt  Fdata(l.l>>  c2  -  aaa 

TvPE*.  lBLI. lBL2  ' . Cl. C2 

IF! I DREC<  4 ;  N£  6)  CO  TO  1113 

IF  ( I  DR EC  <  5 )  LT  15)  CO  TO  1113 

IF  BOTH  ARE  TRUE.  THE  CAINS  ARE  THE  SAME  FOR 
EACH  PLOT 

IF (Cl  GT  C2>  C2-C1 
TYPE*.  Cl.  C2 
IF< C 1  LT  C2)  C1-C2 
TYPE*.  Cl.  C 2 

1 i  13  IF ( IFIRST.  EQ.  1 ) CALL  PLOT!  5.  2  .-3) 

TITLE- 'ATTENUATION  -  ' 

IF < IFIRST.  Ed.  0 )  CALL  PLOT< 1 1.  . 0.  . -3) 

CALL  SYMBOL! 3,  7.  ,  .  1.  XREF! TITLE) ,  90  .  14) 

CALL  NUMBER < 3,  B.  6,  .  2,  FLOAT! IOREC( 16) > ,  90  .  1 > 

CALL  PLOT ( 0.  ,0  .3) 

CALL  PLQTT (FTIM.  FDATA,  3.  XST,  SCALE.  5000.  IREC,  el,  c2> 
IFIRST-0 

C 

C  ADD-UP  NUMBER  OF  DATA  POINTS  EXCEEDING 

C  AMPLITUDE  X  AND  DIVIDE  fiY  TOTAL  DATA 

C  POINTS  IN  RECORD.  USE  LOGARITHMIC 

C  PLOTTING  INTERVAL 

C 

NPT-100 
A  -  5000 
B  -  5000. 

DO  401  ICH-1.2 
C-AVEKK 

IF! ICH  EG  2)  C-AVELL 

DO  351  K-l.npt 

FX-K*C 

X ( K ) -F  X 

COUNT  -  0 

DO  501  L- 1.5000 

IF  <FDATA< ICH.  L ) .  GT.  FX )  COUNT  -  COUNT+1 

501  CONTINUE 

PERCENT ( ICH.  K)  - 
l  (PERCENT (ICH.  K) #A+COUNT ) /B 

WRITE! 6. 1 )  ICH.  K<  X ( K ) .  COUNT,  PERCENT ( ICH, K) 

1  FORMAT ( '  CHANNEL, K,  X.  COUNT, PERCENT  '. 215, 3F15  4) 

351  CONTINUE 

401  CONTINUE 

CO  TO  1005 

900  WRITE( 6,  2) 

2  FORMAT!'  ERROR  IN  INITIALIZING  TAPE') 

GO  TO  1000 

901  WRITE( 6. 3) 

3  FORMAT!'  END  OF  TAPE  ASSUMED') 

PLOT  THE  CUMULATIVE  APD  FOR  THE  ENTIRE  TAPE 
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CH2»C2/80 
XSCAL£»3.  333 
YSCALE-1.  333 
CALL  PLOT < 11  0.  0  , -3) 

CALL  P LOTT <  X .  PERCENT,  1,  XSCALE, YSCALE, npt,  IREC. CHI, CH2) 

CO  TO  1007 

WRITE! 6. 1006)  IREC 

FORMAT ('  END  OF  DISK  OR  TAPE  FILE  DETECTED 
l  AFTER  RECORD  NUMBER  ',15) 

NOW  PLOT  TAPE  ID  PARAMETERS  ON  PLOT  HEADER 

IDREC  <  29) « IREC 
IDREC ( 28)  «  IPAR 
IDREC! 27)  -  1 
CALL  PLOT <9  5,  0  .-3) 

CALL  HEADER (NID) 

CALL  PLOTND 
CLOSE! UNI T-l) 


n  n  c»  o  c\  r j  n  ft  o  o  r<  r»  o  r  i  r>  u  o  o  o  r>  » .  <*j  r»  o  o  r» 


PROGRAM  TPD 


THIS  PROGRAM  COMPUTES  THE  TIME  PROBABILITY 
DISTRIBUTION  FOR  A  USER-SELECTED  SAMPLE  OF  DATA 

WRITTEN  BY  DR  STEN  QDENWALD 


COMMON /HEADER/  IDREC<30>. 1ST (4), STSC, SPSC 
COMMON/DATA/  IDATA<2,  10000) , NCH 
COMMON/AUTO/  TPD (2- 3. 5000) 

COMMON/TIME/  IH.  IM,  SEC,  ATIM< 5000) 

CHAR AC TER *30  FILENAM 
CHAR AC  TER* 14  TITLE 
CHARACTER  BUFFER *32768 
BYTE  J DATA <  32768 ) ,  TEST ( 100 ) 

EQUIVALENCE  < JDATA, BUFFER > 

DIMENSION  FMAX<2, 2000),  *on»< 5000) , REF1 < 100) ,  REF2( 100) 
DIMENSION  ( tpd (2, 3>,  X ( 5000) . PERCENT (2,  5000). SIGMA(2.  2000) 
DIMENSION  AVE <  2,  2000) .  FTIM<5000>.  FDATA12, 5000) , TIM(2000) 
LOGICAL  ERRFLG 
TYPE*.  'TYPE  INCREMENT' 

READ*, INNX 
NX-1000/INNX 
DO  8870  1*1, NX 
X( I )»INNX*< 1-1 ) 

9870  CONTINUE 

NUMBER  OF  STEPS  IN  AUTOCORRELATION  FUNCTION 

NX-32 

NX*10 

OPEN  THE  TAPE  FILE  ON  THE  SCRATCH  DISC  CONTAINING 
THE  DATA  FROM  THE  GWEN  OBSERVING  RUN 

TYPE*,  'INPUT  TAPE  FILE  NAME' 

READ <5, 333)  FILENAM 
F I LENAM* ' S I A 1 :  CODENWALDIB. DAT' 

33  FORMAT < 5X. A30) 

OPEN  <  UN IT- 1 ,  CARR  I ACECONTROL- 'LIST', 

1 S  T A TUS- '  U WNOWN  ',  RECORDS  I ZE-20040) 

INITIALIZE  PLOT  UTILITY 
CALL  PLOTST < 5,  'IN',0> 

IFIRST-1 

IDATA  CONTAINS  THE  INTEGER  VALUES  FOR  UP  TO 
10.000  DATA  SAMPLES  FOR  TWO  CHANNELS  EACH 
HAVING  N1  AND  N2  SAMPLES  RESPECTIVELY. 

ITPD  IS  THE  VECTOR  CONTAINING  THE  THREE 
THRESHOLD  SENSITIVITY  LEVELS  FOR  THE  TPD 


PROGRAM  BEGINS 

DEFINE  NUMBER  OF  WORDS  IN  DATA  RECORD 

ID  RECORD 

NOA  -  l 0000 
NID  -  20 
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NTOT-NDAfNID 
NS  Y  «  2*NT0T 

NUMBER  OF  TAPE  RECORDS 

MREC  »  1120 
NSKIP  ■  0 
DO  88  l-l ,  30 
IDREC  < I )-0 
CONTINUE 
DO  23  1-1.2 
DO  24  J-l, 1000 
PERCENTU.  J)-0. 

CONTINUE 

CONTINUE 

type*.  'record  number 
READ*.  Irec 
NSKIP-IREC-1 
MREC-1 
IDREC( 30) -1 

READ  TAPE  RECORD  OF  LENGTH  NBV  BYTES 

ISTT-1 
IPAR  -  0 
NCH  -  2 
IPLOT-1 

IF (NSKIP.  EQ.  0)  GO  TO  202 
DO  201  IX  — 1.  NSKIP 
READ!  1.5)  BUFFER 
FORMAT (A) 

CONTINUE 

READ (1.3.  END- 1 004 )  BUFFER 
WRIT£( 6. 888)  (JDATA(I),  I-l.  10) 

FORMAT (1015) 

CALL  DECODE <  JDAT A.  NBY,  IREC.  ISTT) 
see  if  the  gain  setting  has  changed 
GAIN-1. 

JCAIN- IDREC (16) 

IF(JCAIN.  EQ.  0)  JCAIN-t 
IF (JCAIN.  EQ.  2)  GAIN-  0631 
IF(  JCAIN.  EQ.  4)  GAIN  -  .00668 
IF( JCAIN.  EQ.  8)  GAIN  -  .000708 

DELT-1000.  76200. 

SUM1-0. 

SUM2-0. 

KK-1 

DO  32  1-1.5000 
FTIM( I )-I*DELT 

fdatad,  i  )-<  l.  0*idata(l.  i)/16.  )**2. 
f data ( 2.  i>«(l.  0*IDATA(2.  I>/16.  )**2. 
SUM1-SUM1+FDATA( l, I) 

SUM2-SUM2+FDAT A ( 2.  I) 

CONTINUE 

AVEKK-SUm/5000. 

AVELL - SUM2 / 5000 . 

XST-10. /( 5000.  *DELT ) 


n  r> 


VSl-  1/AVERK 
VSO  »  1/AVEU- 

IfUnnin*  the  y  ms  limit  IN  UNITS  OF 
10<-13>  WATTS/METER  SQUARE/HERTZ 

SCALE-YS1 

IF(FDATA< 1,  10;  LT  FDATA ( 2. 10) )  SCALE-YSO 
AAA-8  *8.  62*-4/<GAIN*SCALE) 

BBB-S  *8.  62*-4*i0**< 30.  '20.  > / < GAIN*SCALE ) 

cl  ■  im 

c2  *  bbb 

if< #data(2.  1 i  gt.  fdatad.  1  >  >  cl  »  SbB 
iH  fdata(2.  1 )  gt  fdatad,  1>)  c2  »  uj 
TYPE*.  'LBL1.  LBL2'.  Cl.  C2 
IF( IDREC ( 4 ) .  NE  6 )  SO  TO  1113 
IF ( IOREC ( 5  >  LT  15)  GO  TO  1113 

IF  BOTH  ARE  TRUE.  THE  CAINS  ARE  THE  SAME  FOR 
EACH  PLOT 


IF  (Cl.  CT.  C2>  C2-C1 
TYPE*,  Cl,  C2 
IF<C1.  LT.  C2)  C1-C2 
TYPE*,  Cl,  C2 

IF<  IFIRST.  E<3.  1>CALL  PLOT!  5,2  , -3) 

TITLE* 'ATTENUATION  *  ' 

IF( IFIRST.  EQ.  0)  CALL  PLOTd  1.  .  0.  , -3) 

CALL  SYMBOL!-.  3,  7  ,  1,  ‘/.REF (TITLE! .  90.  .  14) 

CALL  NUMBER ( -  3,  8  6,  2,  FLOAT ( JCAIN) ,  90  .1) 

CALL  PLOT  ( 0.  .0.3) 

CALL  PLOTT (FTIM,  FDATA.  3.  XST,  SCALE.  5000.  IREC. cl.  c2) 
IFIRST -0 

TYPE*.  'TYPICAL  SIGNAL  VALUES  ' 

TYPE*.  fdatad,  5).  fdata<2. 5) 

TYPE*,  'TYPE  TWO  THRESHOLDS  TO  USE' 

READ*.  FTPDd.  2).  FTPDd.  3) 

READ*.  FTPD ( 2,  2 ) .  FTPD(2.  3) 

FTPDd.  D-l.  0 
FTPD (2. 1)*1.  0 

PLOT  THE  THRESHOLDS  ON  THE  RECORD  PLOT 

DO  350  1-1.2 

XO-4.  *  I 

DO  351  J-l.3 

XX— XO-SCALE*FTPD( I, J> 

CALL  PLOT (XX.  0  .3) 

CALL  PLOT CXX,  10..  2) 

CALL  NUMBER (XX,  10.  3,  .  1. FLOAT(U>, 90.  ,  1 ) 

CONTINUE 

CONTINUE 

TPINC-1. 

coaputa  th*  d i f f *r tn t ia 1  tpd 

THAT  IS.  THE  NUMBER  OF  TIMES  THAT  THE  PULSES 

ARE  SPACED  BY  I  MILLISECONDS 


'.V.sl 

'  .-’.'vs  m 


DO  700  L-l . 3 
DO  400  J-l.  2 
DO  4,00  1-1.  NX 


•-V-W.- 


77.' 


TP  «  0 

IMAX  =  5000- X i I ) 

DO  550  K-l. I MAX 
INDX-X <  I)+K 

IFlfdata! J. K>  LT  FTPD<  J.  L  >  )  GO  TO  550 
IF!fdata!J. INDX)  LT  FTPD! J. L) )  GO  TO  550 
TP-TP+1 

550  CONTINUE 

IF ( I  EG. 1)  TPO-TP 
TPD<  J.  L,  I  >«  TP/TPO 

;  TYPE*.  'TP,  TPO,  TPD<  J,  L,  I  )  '.  TP,  TPO.  TPD(  J,  L.  I  ) 

ifCtpdl j. 1.  i)  It  .00001)  tpd<j,  1.  1 >-  00001 
c  WRITE! 6. 93)  J.  L,  I .  X < I ) ,  FTPD! J,  L ) , TPD< J. L,  I ) 

:  3  3  FORMAT!'  J,  L,  I ,  X,  TRSH,  TPD  '.315,  3F15  5) 

aOO  CONTINUE 

c  now  sum  the  differential  tpd  to  find  the  integral 

•:  IE  THE  NUMBER  OF  TIMES  THAT  THE  PULSE  SPACINGS 

C  ARE  GREATER  THAN  I 

NXX— NX— 1 

>:  TYPE*.  'CHANNEL,  THRESHOLD  ',  J,  L 

do  601  i-1. nxX 
sum!  i  )  >0.  0 
d o  602  i  i  «  i  ♦  1 . ns 

IF(TPD ! J.  L.  I)  LT.  .  00001)  TPD( J. L,  I)-.  00001 
sum! i ) -sum! i )+tpd ! j. 1. ii) 

1  TYPE*.  'I.  SUM',  I,  SUM!  I) 

IF! I  NE. 1 )  GO  TO  602 

C  COMPUTE  NORMALIZATIONS  TO  INTEGRAL 

C  OVER  ALL  TIME  INTERVALS 

IF!J.  EQ.  1)  REF 1 ! L ) “SUM ! I ) 

IFfJ  EQ. 2)  REF2!L)*SUM! I ) 
c>02  continue 

o03  IF!J  EQ. 1) 

1  TPDI-  SUM! I ) /REF1 I L ) 

IF <  J.  EQ.  2) 

1  TPDI-  SUM!I)/REF2!L) 

C  TYPE*.  'I.  SUM! I).  TPDI'.  I.  SUM! I ) . TPDI 

IF !ABS !TPDI-1.  ) .  LT.  .001)  TPDI-.  977240 
IF! TPD I  CT.  .  977240)  TPDI-.  977240 
IF ! TPDI . LT.  0000001 )  TPDI-1.  OE-10 
TP D!U.  L,  I)-l  -AL0G10!AL0G10! 1.  /TPDI) ) 

C  TYPE*.  J.  L.  I.  SUM!  I).  TPDI.  TPD!  J,  L.  I> 

C  IF! J. EQ.  1 )  WRITE!6.  3334)  I,  J.  L,  SUM! I ) ,  REF1 !L> , TPDI , TPDI 

C  IF!J.  EQ.  2)  WRITE <6. 3334)  I,  U.  L,  SUM! I ) ,  REF2!L) . TPDI, TPDI 

3334  FORMAT!'  I* J,  L,  SUM,  REF,  LOGSM.  TPD  ',  315.  4F15.  5) 

601  continue 

400  CONTINUE 

700  CONTINUE 

GO  TO  1005 

900  WRITE! 6, 2  > 

2  FORMAT!'  ERROR  IN  INITIALIZING  TAPE') 

GO  TO  1000 

901  WRITE! 6, 3) 

3  FORMAT!'  END  OF  TAPE  ASSUMED') 

C 


PLOT  THE  TPD 


:0G5  CHI  *  <  1 0  /5000  )  * l 000  *  8057 

CH2-C2/80 
xSCALE-1,  333 
VSCALE-1.  333 
CALL  PL0T(  11.  0,  0  .  -3) 

WR  ITE  (  6,  3333) .  <  (  ( TPD<  I,  J,  K>.  1*1,2).  J=  1,3),  K-l .  NXX  ) 

3*33  FORMAT (  '  '-6F10  5) 

CALL  PLOTT  (X,  PERCENT,  2,  XSCALE,  VSCALE.  mX,  IREC,  CHI,  CH2) 
GO  TO  1007 

.304  UR  I T£( 6,  1006 i  IREC 

■  ' Oe  FORMAT (  '  END  OF  DISK  OR  TAPE  FILE  DETECTED 
1  AFTER  RECORD  NUMBER  ',15) 

NOW  PLOT  TAPE  ID  PARAMETERS  ON  PLOT  HEADER 

-C0~  IDREC  <  29) » IREC 

IDREC <  28)  «  IPAR 
IDREC (27)  «=  1 
CALL  PLOT (9  5, 0  , -3) 

CALL  HEADER (NID) 

.000  CALL  PLOTND 

CLOSE < UNI T*1 ) 

STOP 

END 


.  TITLE 


BYSWAP  -  SWAP  BYTES  OF  WORDS  AND  MASK  LOWER  10  BITS 


SUBROUTINE  TO  SWAP  BYTES  IN  A  16  BIT  WORD 
INPUT: 

4(AP>  -  ARRAY  ADDRESS  OF  16  BIT  ELEMENTS  TO  BE  BYTE 
SWAPPED  AND  MASKED 


t8(AP) 

OUTPUT: 


NUMBER  OF  16  BIT  ELEMENTS  IN  THE  INPUT  ARRAY 
TO  PROCESS 


LOOP: 


RETURN: 


THE  RESULTANT 

ARRAY  OF  16  BIT  SWAPPED 

ELEMENTS 

*  •  -  -  a 

OVERLAYS  THE 

INPUT  ARRAY 

w.  , 

SUBROUTINE  IS 

FORTRAN  CALLABLE.  AS  IN: 

CALL 

BYSWAP ( I ARRAY.  NUMBER ) 

yy-y/y- 

WHERE  I ARRAY 

IS  THE  1*2  ARRAY  OF  DATA 

TO  BE  PROCESSED  AND 

NUMBER  IS  THE 

1*2  NUMBER  OF  16  BIT  ELEMENTS  IN  I ARRAY. 

.9, 

. PSECT 

CODE.  RO.  NOWRT.  EXE.  LONS 

.  ENTRY 

BYSWAP.  •'M<IV.R2> 

MOVL 

4(AP>.  R1 

] PICK  UP  ARRAY  ADDRESS 

CLRL 

R2 

MOVW 

•8 ( AP ) .  R2 

(R2  -  NUMBER  OF  ELEMENTS  TO  PROCESS 

BLEQ 

RETURN 

t DO  NOTHINC  IF  INVALID  COUNT 

• 

EXT7V 

*9. *8.  (Rl).RO 

»  PLACE  8  MSBits  of  iourci  into  LSbits  of  RO 

INSV 

<R1>.«8.  «8.  RO 

land  8  LSBitt  of  lourct  into  MSbit»  of  RO 

MOVW 

RO.  <R1  )■*■ 

imov*  mult  to  ovirwrit*  input 

-•%V. 

SOBOTR 

R2. LOOP 

i do  again  for  **ch  *l»m»nt 

>  •?  ’* 

RET 
.  END 

(DONE,  RETURN  TO  CALLER 
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